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A canonical formalism and constraint analysis for discrete systems subject to a varia- 
CN| ■ tional action principle are devised. The formalism is equivalent to the covariant formulation, 

encompasses global and local discrete time evolution moves and naturally incorporates both 

■ constant and evolving phase spaces, the latter of which is necessary for a time varying dis- 
ij cretization. The different roles of constraints in the discrete and the conditions under which 

they are first or second class and/or symmetry generators are clarified. The (non-) preser- 
vation of constraints and the symplcctic structure is discussed; on evolving phase spaces the 
number of constraints at a fixed time step depends on the initial and final time step of evolu- 
tion. Moreover, the definition of observables and a reduced phase space is provided; again, on 
evolving phase spaces the notion of an observable as a propagating degree of freedom requires 
£SJ ' specification of an initial and final step and crucially depends on this choice, in contrast to 

^ . the continuum. However, upon restriction to translation invariant systems, one regains the 

usual time step independence of canonical concepts. These results are applicable, e.g., to 
. discrete mechanics, lattice field theory, quantum gravity models and numerical analysis. 

■ 1 Introduction 

O 

m 

Discretizations of dynamical systems assume a pivotal role in both numerical simulations of 
physical models and the non-perturbative quantization of field theories. Quantum gravity is a 
prominent example of the latter in which discrete structures appear in many different ways. 

Usually, it is space—time itself which is discretized, for instance, by a triangulation. In gravity 
space-time is a dynamical entity, such that one may expect some aspects of this discretization 
also to be dynamical, i.e. to change in time. For instance, different hypersurfaces in a generic 
space-time triangulation in Regge Calculus [1] — the most well-known discretization of General 
Relativity — are comprised of different numbers of simplices [2]. In particular, one may wish to 
model an expanding/contracting universe with a time dependent discretization, in which the 
number of degrees of freedom grows/shrinks with the size of the spatial geometry O d]. In 
the quantum theory, this leads to the notion of 'evolving Hilbert spaces' [5]. In the numerical 
analysis of physical systems, on the other hand, it is often convenient or (for computational 
savings) even necessary to adapt either to a more refined or coarser discretization during the 
evolution of the system. This, in particular, has culminated in the standard tool of 'adaptive 
mesh refinement' [6]. 
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In such situations one has to cope with discrete systems whose numbers of degrees of freedom 
change in time. In order to systematically treat such scenarios, first of all at the classical level, 
the present work shall develop a general (classical) canonical formalism applicable to arbitrary 
discrete systems governed by a variational action principle. In particular, this formalism applies 
to discrete mechanics, lattice field theories and certain discrete gravity models. Its general 
features are the following: 

• It can handle both constant and evolving phase spaces, the latter by natural phase space 
extensions. 

• It is equivalent to the covariant formalism. By using the action (or Hamilton's principal 
function) as a generating function, the canonical formalism is directly derived from the 
covariant formulation. Furthermore, ensuring that all constraints (on the extended phase 
spaces) are always satisfied is tantamount to solving the (covariant) equations of motion. 

• It is insensitive to the particular discretization and form of the (effective) action. Hence, 
it is amenable to coarse graining methods (which are relevant for the continuum limit). 

We shall see that systems with changing phase space dimension necessarily lead to con- 
straints. It is essential to distinguish the latter type of constraints from symmetry generators, 
in particular, in theories of gravity. In continuum (Einstein) gravity, the diffeomorphism sym- 
metry leads to the Hamiltonian and spatial diffeomorphism constraints, which also generate 
(infinitesimal) time evolution. The situation in the discrete is, however, very different because 
diffeomorphism symmetry is generically broken [TJ El El [10] such that Hamiltonian and diffeomor- 
phism constraints (generically) do not arise. Rather, equivalence of the canonical and covariant 
picture requires canonical time evolution to proceed in discrete steps [2]. Nonetheless, a changing 
phase space dimension will always lead to constraints, which only in special situations include 
or even coincide with the Hamiltonian and diffeomorphism constraints. 

A proper understanding of the role of the constraints is thus important for the interpretation 
of both the classical and quantum theories based on such discretizations (SHU [11]. In this article, 
we shall provide the necessary constraint analysis for variational discrete systems which in some 
aspects is analogous to and in others quite dissimilar from the continuum Dirac procedure 
[12^13] . In the continuum, the roles of the constraints pertain, in particular, to the identification 
of gauge symmetries and propagating degrees of freedom, as well as observables. On evolving 
phase spaces, one has to reconsider the very concepts of propagating degree of freedom and 
observable; as we shall see, their definition, in contrast to the continuum, cannot be associated 
to one particular instant of time only. 

The constraint analysis for variational discrete systems has the following main results: 

• The (non-) preservation of constraints and the symplectic structure on evolving (and 
constant) phase spaces is clarified. In general, the number of constraints at a fixed time 
step critically depends on the initial and final steps of evolution because equations of 
motion can lead to a 'propagation' of constraints. 

• The conditions under which the canonical constraints are symmetry generators and first 
or second class are investigated. 

• The concept of observables as gauge invariant and propagating degrees of freedom requires 
two time steps on an evolving phase space and crucially depends on these initial and final 
steps under consideration. For two different pairs of steps one, in general, finds a different 
number of propagating degrees of freedom. 
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• The reduced phase space at some step n in an evolution i — > n — > / depends on initial and 
final steps % and /. The classification of degrees of freedom is, in general, step dependent. 

The step dependence of many canonical concepts notwithstanding^] the new formalism is fully 
consistent, solves the covariant equations of motion and unambiguously describes the propaga- 
tion of data between different time steps. It is therefore applicable to evolving lattices and offers 
a comprehensive picture of the discrete dynamics. 

For example, as shall be explained in the course of this article, systems which have been 
evolved from a phase space with zero degrees of freedom (i.e. from 'nothing' as in the 'no 
boundary proposal' [13]) are always totally constrained, in the sense that the dimension of the 
reduced phase space on the evolving 'time slice' is zero. On account of the step dependence of 
observables, however, this does not mean that these systems are necessarily devoid of propagating 
degrees of freedom. 

These results may prove useful for numerical implementations in discrete mechanics and 
lattice field theory and shed light on the interpretation of quantum transition amplitudes, e.g., 
in loop quantum gravity [15|. [T6l \17\ . Moreover, they further the interpretation of the newly 
developed canonical formulation of Regge Calculus [21 [18] . These frameworks also involve dis- 
cretizations (graphs), encoding different numbers of degrees of freedom for the initial and final 
state. In fact, the present formalism shares certain features with the 'general boundary formu- 
lation' \19\ I20|. where transition amplitudes are defined between very different types of surfaces 
(including 'zero-size' surfaces) and the traditional splitting of the space-time manifold into 
A4 = S x R is not assumed. 

The rest of this manuscript is organised as follows: in section [2[ we shall briefly introduce 
variational discrete systems and, subsequently, discuss the global dynamics of such systems in 
the regular and irregular case in section [3] (for a detailed exposition of the regular case see 
[21)). The irregular setting encompasses systems with evolving numbers of degrees of freedom. 
In particular, section [3] introduces the Legendre transforms and the canonical formalism for 
discrete systems. We shall provide sufficiently many details to make the present article relatively 
self-contained. 

Section 2] concerns local evolution moves in which only a (small) subset of an equal time 
hypersurface is evolved. A typical example are the Pachner moves which constitute the most 
elementary evolution steps for a triangulated equal time hypersurface. For Regge Calculus, these 
have earlier been introduced as canonical evolution moves in [2], below we shall generalize this 
concept to general variational discrete systems. Section [5] analyzes the different roles constraints 
can assume in the discrete and specifies in which cases they constitute gauge symmetry genera- 
tors. Moreover, we generalize the concept of observables and propagating degrees of freedom to 
discrete systems with varying numbers of degrees of freedom. We shall identify such observables 
and provide a counting of the number of propagating degrees of freedom. Finally, we shall close 
with a summary and discussion in section (6[ 

Along the way, we shall illustrate general features by means of a scalar field living on a 2D 
triangulation. Some technical details have been moved to the appendices. 

2 Preliminaries: variational discrete systems 

In this manuscript we shall consider general variational discrete systems whose evolution pro- 
ceeds in discrete time steps which we label by n € Z. 

1 Upon restriction to translation invariant systems this step dependence disappears. 
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Let Q n be the configuration space of the system at time step n. Q n shall be coordinatized by 
x l n , where i takes value in some index set determined by the dimension of Q n . (We shall often 
omit the index i for notational convenience.) We expressly emphasize that the Q n at the various 
steps n need not be of the same dimension. The dynamics of such systems shall be described 
by discrete actions of the form 

K 

SK = y^ J S n {x n -\ , x n ) , (2.1) 

n=l 

where the sum ranges over the individual time steps n. The individual action contribution S n 
governs the discrete time evolution move (n — 1) — > n, i.e. the discrete time evolution from time 
step n — 1 to time step n. These systems are 

(i) variational because the configuration spaces Q n are continuous manifolds such that the 
dynamics is gained from a variational action principle, 

(ii) discrete because the time evolution proceeds in discrete steps labeled by n, and, finally, 

(iii) mechanical because Q n shall be finite dimensional. 

The main assumption we shall henceforth make is that of additivity of the discrete action, 
in order for the sum (|2.ip to make sense. That is, the action contribution of the union of any 
two evolution moves should simply be the sum of the individual action contributions of the two 
moves. In particular, if the two moves {n — 1) — >• n and n — > (n + 1) share a 'common boundary' 
described by the variables x n , we can solve the equations of motion 

dS n (x n —i,x n ) dS n +i(x n i x n +\) 
dx n dx n 

for x n = Xn{xn-lt x n+i) i n order to obtain a new 'effective' action contribution 

S(n—l)->(n+l)[ x n—li x n+l) := Sn \ x n— li Xn\ x n— It x n+l 

)) + S n+ i (x n ( x n-l, x n+l), x n+l) , (2-3) 

governing the 'effective evolution move' (n — 1) — > (n + 1). The action evaluated on solutions, 
given suitable boundary data, is called Hamilton's principal function (2TJ [22], [2j [231 EH [25] and 
is a function of this boundary data. In our case Hamilton's principal function Sr n _i\_^r n+ i\ is a 
function of £ n _i and x n+ i coordinatizing the initial and final boundary, respectively, and defines 
a new contribution in the action sum (|2.ip that summarizes two elementary ones. That is, for 
full generality we allow for the two possibilities that a given individual contribution S n in (I2.ip 
may either be an elementary ('bare') action or Hamilton's principal function that summarized 
several elementary evolution moves into a single effective one. 

Such variational systems encompass a large variety of discrete physical systems, ranging from 
discrete mechanics [21], over lattice field theory [26] to discrete gravity. In particular, Regge 
Calculus falls into this class of discrete mechanical systems [2] [18]: the lengths of the edges of a 
Regge triangulation can take continuous values, in a finite triangulation there are finitely many 
edges and the discrete time evolution proceeds by discrete evolution moves (e.g. Pachner or tent 
moves [2]). Furthermore, due to boundary terms, the Regge action is additive in the above sense 

[23 E] IT 

2 The Regge action is given by Si^ggc = J2 hc7 -\ gj - VuSh + J^hcdT^^^ 1 ^ wnere h denotes the hinges, i.e. 
(D — 2)-dimensional subsimplices in the D-dimensional triangulation T, Vu denotes the volume of h, Eh denotes 
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(2.2) 



Most research on such systems has been performed in the Lagrangian setting. It is the 
goal of the present work to develop a canonical formalism for such general variational discrete 
systems which is equivalent to the Lagrangian formulation. This is achieved by using the action 
(or Hamilton's principal function) S n = S n (x n -i, x n ), which depends on 'old' variables x n -\ 
and 'new' variables generating function of the first kind for canonical time evolution 

|2H [22| [23| [2^1 [21 [HI [18] . The new canonical formalism shall be applicable to general systems with 
evolving phase spaces and, in particular, provide a constraint analysis for variational discrete 
systems in analogy to Dirac's original continuum procedure |12|, |13| . Finally, given that the 
formalism naturally handles Hamilton's principal functions, different choices for the underlying 
discrete lattice can be related to each other and the formalism is naturally applicable to coarse 
graining methods (which, e.g., are relevant for the continuum limit), see for instance [28 ^ I29| , [30] . 

Example 2.1. An example which we will use throughout this paper to illustrate the general 
principles will be the 2D discretized massless scalar field on an triangulation with Euclidean 
geometry. To simplify the action we shall consider an equilateral triangulation. The scalar field 
<p will be associated to the vertices of the lattice, thus <f> v denotes the field at the vertex v. 
The action associated to one equilateral triangle is given by |31] 

sa = \ ( £(<n 2 - E ^ s(f V (e) ) = i E(^ (e) - ^ t(e) ) 2 ' ( 2 - 4 ) 

ViiCA ecA / ecA 

where v, e denote the vertices and edges of the triangle A and s(e), t(e) are the source and target 
vertex of the edge e. The action associated to any larger triangulation is given by the sum of 
the actions associated to the triangles. This automatically includes the right boundary terms, 
if the fields on the boundary are held fixed. 

3 Global dynamics of variational discrete systems 

In this section we shall only consider global time evolution moves. Global moves are such that 
each of the variables at a given discrete time step is involved in the move and only occurs at this 
one time step, i.e. neighbouring time steps n, n + 1 do not overlap, except for variables in the 
boundary. For example, in simplicial gravity this corresponds to evolution moves which evolve 
between disjoint spatial hypersurfaces [2]. Local evolution moves which only evolve subsets of 
data of a given time step will be discussed in the subsequent section UJ 

3.1 General properties and definitions 

Consider three consecutive steps n — 1, n, n + 1 and the boundary value problem defined by the 
data at times n — 1 and n + 1. That is, we are given boundary data x t n _ 1 and x l n+1 and ought 
to extremize 

S '■= S n (x n —i, x n ) + S' n _)_i(x n , x n +i) (3-1) 

the deficit angle around hinge h in the bulk T \ dT of the triangulation and ?/>h denotes the extrinsic curvature 
angle around hinge h in the boundary dT of the triangulation. When two pieces of triangulation Ti and Ti are 
glued together, (parts of) their boundaries and, in particular, the hinges in the boundaries are identified. The 
two extrinsic curvature angles around a fixed hinge h in both boundaries of 71 and Ti add up to a deficit angle, 
V'hceTi + V , hcer 2 = £hc(7iuT 2 )\8(7iuT9)! i n tne bulk of T\ U 72- The Regge action is therefore additive. On the 
other hand, a discrete action where Sh is replaced by sin Eh constitutes an example of a non-additive action. 
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with respect to x n . This yields the equations of motion 

= 9Sn | dSn+1 t (3.2) 

which may or may not be uniquely solvable for x n as a function of x n —i,x n +i, depending on 
whether the system under consideration is regular or irregular (see sections 13.21 and 13.31 below) . 



3.1.1 Lagrangian formulation 

The action contribution S n defines a mapping S n : Q n -i x Qn —> On the other hand, 
in continuum mechanics, the Lagrangian L(q, q) defines a mapping L : TQ —¥ K. That is, 
in the discrete the direct product of configuration manifolds Q n -i X Q n — coordinatized by 
x n _i,x n — assumes the role of the tangent bundle TQ — coordinatized by q, q — of the Lagrangian 
formulation of the continuum. In fact, if Q n -i — Qn — Q, Q x Q is locally isomorphic to TQ. 

The variation of the discrete action (13.1j) enables one to define the so-called Lagrangian 
one- and two-forms [21]. In contrast to the continuum where only one Lagrangian one-form 9 
exists on TQ, in the discrete two Lagrange one-forms on Q n x Q n +i arise from the boundary 
terms of the variation of the action. Namely, varying S n +\ as 5S n+ \ = dS n +\ ■ S q . +1 . with 
q n ( n +\) G Qn x Qn+i and some variation 5 Qn{n+1) £ T qn{n+1) (Q n x Q n +i) yields 

dS n+ i = n+1 - 9~ , (3.3) 
where (summing over repeated indices i,j is understood) 

a—/ \ _ d&n+l j i 

v n \ x rii x n+l) — a j a X n 



dx l n 



9+ +1 (x n ,x n+ i) = 9S ™ +1 dx n+1 . (3.4) 

OX n+l 



However, d o dS n+ \ = 0, such that a single Lagrange two-form on Q n x Q n +i can be defined: 



d 2 S n 

n n+ i(x n ,x n+1 ) = -d6+ +l = -dB~ = - . "+ 1 dx\ /\dx J n+1 . (3.5) 

oxlox J n+1 



The equations of motion (|3.2p can be used to define a Lagrangian time evolution map C n '■ 

Qn-l X Q n -)■ Q n X Q n+1 

£-n '■ ( x n—l,%n) l— ^ (%n> x n+l)> (3-6) 

by solving for x n +i given (x n -i,x n )- C n may neither be defined on all of Q n -\ x Q n , nor map 
to all of Q n x Qn+i, nor be unique in the presence of constraints. This shall be a subject of the 
following subsections 13.21 and [ 



3.1.2 Discrete Legendre transformations and canonical formulation 

In order to discuss the dynamics in a canonical language, we need to introduce discrete Legendre 
transformations which will carry us to suitable phase spaces. Recall that in the continuum 
formulation a single Legendre transformation FL : TQ — > T*Q exists, which in coordinates 
reads (q,q) h4 (q, On the other hand, two Legendre transformations can be defined in the 
discrete because we work with a cartesian product of two configuration manifolds Q n x Q n +i, 
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instead of TQ. The precise definition of these discrete Legendre transformations which we will 
employ is motivated by, however, suitably differs from the continuum version (for details on the 
continuum Legendre transform see |22j). 

Consider an arbitrary step n. Recall that S n : Q n -i x Qn ^ where Q n -\ x Qn is a 
fibre bundle. Pick a point q n -i € Q n -i (see figured]). We denote the fibre over fixed q n -i by 
F n (q n -i) := (Qn-i x Qn)q n -i- Notice that T n — Qn- Choose a point f n € T n and a curve 7(e) 
in T n with curve parameter e such that 7(0) = /„ and 7 7 (0) = ^7(e) |e=o- This allows us to 
provide the following 

Definition 3.1. TTie discrete fibre derivative ¥ + S n : Q ra _i x Q n — > T*Q n , defined by 



F + 5„(/„)-7'(0) :=^ n ( 7 ( £ )) 



£=0 



(3.7) 



is called the post-Legendre transform. 

F + 5 n (/ n ) • 7'(0) is the derivative of S n along the fibre T n at / n in the direction 7 ; (0). Given 
that T n = Q„, we have 7 '(0) £ T fn Q n and, thus, F+5„(/„) G r; n Q n . 




-^n(gn-l) = Qn 



/n = 7(0) 



Qn-1 X Q 7 



q n -i 



Q 



n-1 



Figure 1: Illustration of the fibre derivative used in the definition of the post-Legendre transform. 



Now exchange the roles of Q n -i and Q n and choose / n -i € Fn-xifln)- Let r/(zv) be a curve 
in -Fn-i such that 77(0) = f n -i- In complete analogy, we give 

Definition 3.2. The discrete fibre derivative ¥ + S n '■ Q n -\ x Q n —> T*Q n „i, defined by 



u=0 



is called the pre-Legendre transform. 

The cotangent bundles V n -i '■= T*Q n _i and V n '■= T*Q n are the phase spaces which we 
will henceforth work with. 

For the remainder of this article, we will need the coordinate form of the pre- and post- 
Legendre transform. Consider the post-Legendre transforms. In a small neighbourhood of f n 
we can write 7(e) = / n + £7'(0). Let (x n -i,x n ) be the coordinates of f n in Q n -i X Q n . Thus, in 
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coordinates 7(e) reads {x n _\,x n + e5x n ) with some Sx n . Inserting this in S n (x n _i, x n ), yields 
(13.71) in the form 



d 



de 



dS n 



—S n (x n -i,x n + e5x n ) =- — 5x n . (3.8) 



e=0 dx n 



Hence, we obtain as in [2T] (note that x n are the coordinates of f n in Q n ), 

F + S n (x n _ l ,x n ) = (^n,^^) • (3.9) 

The coordinate expression for the pre-Legendre transform is derived in complete analogy. In 
general, we write 

¥+S n : s„) ' y (x n +p n ) = [Xn,^A (3.10) 

¥~S n : (x n -i,x n ) i-> (x n _i,"_p n_1 ) = fx n -i,- ^ n ^ \ . (3.11) 

We shall refer to (|3.10p as the post-Legendre transformation and to (|3.1ip as the pre-Legendre 
transformation. 

The reason for the minus sign in the definition of in (|3.11j) is the following: using the 

coordinate forms (|3.10l 13. lip , it is straightforward to check that — in analogy to the continuum — 
the Lagrangian one- and two-forms, (13, 4h and (|3.5p . arise from pulling back the canonical one- 
and two-forms, 9 n = pfdx l n and uo n = dx^Adpf, respectively, with the Legendre transformation. 
That is@ 

et = (¥ + s n ye n , n n = (w + s n yco n , 

9- = (¥-s n+l ye n , n n+1 = (¥-s n+1 )*u> n . (3.12) 

In order to define the Legendre transforms, we employed the discrete action, or Hamilton's 
principal function. By noting that Hamilton's principal function is a generating function of 
the first kind (i.e. depends on the old and new configuration coordinates), one can define the 
canonical discrete time evolution: 

— n— 1 9S n (x n —i, X n ) + n dS n (x n —i, X n ) 

V '■— ^ ) P '■— ^ • {6.10) 

ox n -x ox n 

We shall refer to the momenta ~p ([3. lip as pre-momenta and to the momenta + p (|3.1Up as 
post-momenta. 

Equation (|3.13p defines an implicit global Hamiltonian time evolution map H n -\ ■ V n —\ — > 

U n ^ : (x B _i, "p"- 1 ) -> (x n , V). (3-14) 

Namely, given (x n -i,~p n ~ 1 ), one can use the equation for the pre-momenta in (I3.13p . in order 
to determine x n and, using this result and the post-momenta equation in (|3.13p . one determines 
+ p n . It depends on the presence of constraints whether one can uniquely determine x n or not. 
This will be the subject of the following subsections. 



3 Clearly, the definition is coordinate independent and so, if a new coordinate system x' n (x n ) on Q n is chosen, 
one finds dx' n A dp' n — dx n A dp n . 
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We emphasize that the definition of ~H n — in contrast to the one for C n — itself does not 
require the equations of motion. Nonetheless, the equations of motion are also necessary in the 
canonical picture for the following reason: In analogy to (|3.13p . we could equally well use S n+ i 
as a generating function. Accordingly, at every time step we have both pre- and post-momenta 
(see figure [2]): 

dS n+ i 



P 



dx n 



P 



8S n 

dx n 



(3.15) 



The equations of motion + dSn+1 



for the vari- 



dx n ' dx r , 

ables equivalent to a momentum matching of pre- 

and post-momenta, 



P 



V 



(3.16) 




+ 1 

S n +1 



Sri 



Figure 2: The implicit global discrete 
Hamiltonian time evolution map and pre- 
and post-momenta. 



such that (on-shell) there are unique momenta for the 
variables at step n. Henceforth, we will often omit the 
superindices + and — at the momenta, implicitly as- 
suming that momentum matching holds. 

In summary, we find that the off-shell (or kinemat- 
ical) quantities are defined by the action associated to 
one time step, while equations of motions are imposed 
by momentum matching. 

Example 3.1. Consider the evolution of the massless 
scalar field with action (12. 4p on a two-dimensional cylindrical space time S 1 x [0, 1] , where we 
assign the time n to the lower boundary S 1 x {0} and time (n + 1) to the upper boundary 
S 1 x {1}. We assume that this space time is triangulated such that we have Q n vertices at time 
n and Q n +i vertices at time (n + 1) and such that we have only one layer of triangles as depicted 
in figure El This corresponds to a triangulation as used in the Causal Dynamical Triangulation 
(CDT) programme [32] ; a vertex i at time n can be connected to a vertex j at time (n + 1) by 
one edge. For very 'small' triangulations a vertex i may also be connected to j by 2 edges. This 
defines an adjacency matrix A^ 1 with A™^ 1 = E if vertex i is connected with vertex j by E 
edges, where E = 0,1 or 2. The canonical time evolution from time n to time (n + 1) can then 
be written as 



7T 



n+1 




+ 



'n+1 



2^n+l 2T , n+l' 



(3.17) 



where 7r" are the momenta conjugate to the fields (j) l n and we assume a cyclic labeling of the 
vertices. For instance, if i = Q n then i + 1 = 1 for objects cft^ and 7T l n . 



3.2 Review: Regular discrete systems 

We begin by reviewing regular discrete systems in which — by definition — constraints do not 
arise. The configuration spaces Q n = Q are of equal dimension at every time step n and the 
condition 

d 2 S n+ i 



det 



dx„dx 



n+1 



^0 



(3.18) 
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holds. A detailed exposition of regular discrete systems is given in [21] and so we shall be brief 
on this subject. 



3.2.1 Lagrangian formulation 

In this regular case, given (x n -i,x n ), (|3.2p can be uniquely solved for (x n ,x n+ i) and the La- 
grangian time evolution map C n in (|3.6p is unique. 

The Lagrangian two-form is preserved under the time evolution map 



n n = ow, (3.19) 

where C* n represents the pull-back of C n . This is proven as follows: by substituting for x n+ \ the 
solutions x n+ i(x n -i, x n ) of the time evolution map (|3.6p . the action S in (|3,ip can be written 
as a function on Q n _i x Q n . The exterior derivative of S on Q n _i x Q n yields only boundary 
terms on account of the equations of motion of the inner variable x nj 

,r,/ \ 9S n j dS n+ i ( dx +1 . dx 3 n+l A 

dS{xn-i,x n ) = —■ — dx n _ x + — - — — - — dx n _ l + ——dx n \ 
9<-i dx 3 n+1 \dx\_-L dx\ J 

= -e~_ 1 + C* n e++x- (3-20) 

d o d = and commutativity of exterior derivatives and pull-backs then imply (|3.19p . It is 
straightforward to generalize this argument to an arbitrary time step difference (m,ri2). 



3.2.2 Canonical formulation 

As a consequence of condition (I3.18p . the Legendre transforms (13.101 13. lip are invertible and 
the Hamiltonian time evolution map ~H n -i (|3.14|) generally possesses a unique solution. One can 
check [21] that T-L n -i = F + o C n -\ o (F + ) _1 , where C n -\ is the Lagrangian time evolution map 
(|3.6|) . This, together with the preservation of the Lagrangian two-form (|3.19p can be used to 
show that the Hamiltonian time evolution is symplectic, i.e. preserves the canonical two-form 

w n _i = (Hn-i)*u n . (3.21) 

The situation for regular discrete systems can be diagrammatically summarized as follows: 



Qn-l X Qn — > Qn X Q n+1 
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3.3 Irregular discrete systems 



Next, we examine the general scenario in which regularity condition (|3.18p is violated. Nonethe- 
less, let Qn-i — Qn Vn, such that an equal number of left and right null vectors L l n and R l n+ i, 
respectively, of the Lagrange two-form 

V n &2Sn+1 =0, d2Sn+l Ri +1 =0, (3.22) 

dxldx> n+l dx\dxP n+1 

arises in an open neighborhood in Q n x Q ra +i. (In this section only, a further index labeling 
the null vectors will be ignored for the sake of better legibility.) As we shall see in subsection 
13.41 below, it is possible to formulate systems with temporally varying numbers of degrees of 
freedom in this irregular fashion. 



3.3.1 The Lagrangian formulation 

A specific solution X ri j f -\ clS cl function of (x n _i,x n ) is not uniquely determined because x n+ \ + 
£R n +i then constitutes another (infinitesimally displaced) solution of (|3.2|) . In consequence, 
arbitrariness in the form of free parameters arises. That is, there exists and N parameter family 
of Lagrangian time evolution maps , I = 1,...,N, if there are iV independent null vectors 
Rn+i- A given C n , obtained after fixing N a priori free parameters maps at most onto a 
(2Q — A) -dimensional space where Q is the dimension of Q. Therefore, C n is either not injective 
and/or only defined on some constraint submanifold of Qn—i X Qw 

On account of (I3.22|) . the Lagrangian two-form (13. 5j) is, obviously, degenerate and possesses 
2 A null directions L l n and R l n+ i- Nevertheless, if there are no constraints on <2 n _i x Q n and, 
hence, a given can De defined on all of Q n -\ x Qn, the Lagrangian two-form is preserved 
under the discrete time evolution. The corresponding argument is identical to the one spelled 
out below ([319]) . 

The same holds for the time reverse of C n @ if some null vectors L n _\ are present. It is as 
this stage useful to distinguish between two types of free data: 

Definition 3.3. We refer to a free parameter 

• A n +i at step n + 1 which cannot be predicted by C n , given the Lagrangian data at steps 
n — l,n, as a priori free, and 

• /U n _i at step n—1 which cannot be postdicted by the time reverse of C n , given the Lagrangian 
data at steps n,n + 1, as a posteriori free. 

As we shall S6G in section [5] below, these free parameters Xnj 

[i n and null vectors L n ,R n at 

some n do not necessarily correspond to gauge freedom. 

3.3.2 The canonical formulation 

Recall that in the continuum the Legendre transform ¥L fails to be an isomorphism if and only 
if det ^ Q^igqj ) = 0, i.e. if and only if the Lagrangian two-form is degenerate. In this case a 
single primary constraint surface arises |13j . 

The analogous state of affairs holds true in the discrete. Definitions 13.11 and 13.21 of the 
pre- and post-Legendre transforms are directly applicable to singular systems. In this case the 

4 More precisely, by the time reverse of C„ we mean the map defined by the equations of motion (|3.2[) . given 
suitable 'final data' on Q n x Q n +i- 
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rank of both Legendre transformations is 2Q — N, as a consequence of the N left and N right 
null vectors ()3.22p . The Legendre transformations thus fail to be onto J! and their images form 
(2Q — 7V)-dimensional submanifolds in the two phase spaces V n and V n +i. We shall denote 
these submanifolds by C~ and C^ +1 , respectively. Notice that dimC~ = dimCi. 1 . 

Definition 3.4. The image of the pre-Legendre transform, C~ := Im(F _ SVt+i) C V n , is called 
the pre-constraint surface. The image of the post-Legendre transform, C^ +1 := Im(F + S' n +i) C 
■p n _l_i, is called the post-constraint surface. 

In contrast to the continuum, we obtain two kinds of constraint surfaces in the discrete since 
there are two kinds of Legendre transforms. The pre-constraints governing C~ are satisfied by 
the pre-momenta ~p n by construction because these arise from the definition of the Legendre 
transform. However, the pre-constraints impose non-trivial conditions on the post-momenta 
+ p n at n, resulting from the previous time evolution move (n — 1) — > n. 

If there are N independent left null vectors of the Lagrangian form (|3.22p (in an open 
neighbourhood), the pre-constraint hypersurface will be of co-dimension N in V n and we choose 
N constraints "C™, I = L...,iV, to describe this pre-constraint hypersurface. Thus, this 
constraint set is irreduciblqj [13] . 

The N (irreducible) post-constraints + C" +1 , r = 1,... ,N, defining C^ +1 are automatically 
satisfied by the momenta +p n+1 (for all initial values) after having performed the evolution move 
n — > (n + 1). By momentum matching, these post-constraints will provide non-trivial conditions 
for the pre-momenta ~p n+l of the next evolution move (n + 1) — > (n + 2). 

The pre- and post-constraint surfaces at fixed step n generally do not coincide, C + ^ C~ . 
Making sure that both pre- and post-constraints are satisfied at each step such that we restrict 
to C + n C~ is the non-trivial challenge and will be further discussed in section 15.11 

For singular systems, the Hamiltonian time evolution T-L n (|3.14p can only be defined as a map 
from C~ to C^ +1 . That is, 1~L n : C~ — > C^ +1 . Nevertheless, TL n is generated by the discrete action 
and implicitly defined by (|3.13p . However, as a consequence of the right null vectors in (I3.22p . we 
are no longer able to solve the pre-momentum equation in (|3.13p uniquely for x n +i: therefore, 
as in the Lagrangian picture, we must specify N a priori free parameters AJ^ +1 , r = 1,. . . ,N 
to uniquely determine x n _|_i(x n , X n +i, ~P n )- Moreover, the post-momenta + p n+1 will in general 
depend on the parameters A n+ i. 

Likewise, on account of the left null vectors L n , we can no longer uniquely express x n as a 
function of x n +i, + p n+l only and iV a posteriori free parameters fi l n , I = 1,... ,N, that cannot 
be postdicted via (the time reverse of) H n and the canonical data at n + 1 need be specified in 
order to determine x n (x n+ i, fi n , + p n+1 ). 

Via the momentum matching (|3.16p . the a priori and a posteriori free parameters A n ,u n at 
step n of the Lagrangian and Hamiltonian time evolution maps coincide. We have one a priori 
free A n per post-constraint and one a posteriori free u n per pre-constraint. 

In analogy to the continuum situation, arbitrariness in the form of free parameters A, \x thus 
appears in the canonical discrete evolution of singular systems. But as in the continuum the 
preservation of the constraints may lead to a fixing of some of these parameters A, u. This will 
be discussed in sections I5.1H5.41 

In consequence of the Hamiltonian time evolution 7i n : C~ — > C^ l+l restricting to the con- 
straint hyper surf aces, H n cannot be a symplectic map. However, it defines a pre-symplectic 

5 In analogy to the continuum, F S„+i simultaneously fail to be isomorphisms if and only if the Lagrangian 
two-form H3.5|) is degenerate. 

6 The gradients of the constraints are linearly independent. 
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map: using the embeddings t n : C n — > V n and : C„ +1 — > "P n +i, the canonical two-forms 

w n = da& A dp!?, w n+ i = d< +1 A dpf" 1 . (3.23) 

can be pulled back to two-forms on the constraint surfaces C~ and C„ +1 , respectively. The 
resulting two-forms are pre-symplectic formal and Hamiltonian time evolution preserves these 
pre-symplectic forms as stated in the following theorem. 

Theorem 3.1. The discrete global Hamiltonian time evolution map T-L n : C~ — > C^ +1 satisfies 

(i-)*u; n = n* n (i+ +i ru; n+1 . (3.24) 
Proof. See Theorem 6.1 in [2]. □ 
Let us conclude the treatment of singular discrete systems with a diagrammatic symmary: 



Qn-l X Qn — > Qn X Q n+ i 




n+1- 



3.4 Evolving phase spaces and configuration and phase space extensions 

In order to describe a growing/shrinking lattice, we must be able to deal with numbers of 
variables which vary in discrete time, such that the concept of an evolving phase space becomes 
necessary. 

Note that the Definitions 13.11 and 13.21 of the pre- and post-Legendre transforms are also 
applicable to a situation where Q n -i ¥ Qn- In this case the transforms obviously fail to be 
isomorphisms and we are necessarily in the situation of a singular system. The mapping from 
Qn-l x Qn to the phase spaces V n -i = T*Q n _\ and V n = T*Q n is well defined also in the case 
dimQ n _i ^ dimQ n . Therefore, all the results of section [3TT1 apply to a Q n which varies with 
n. The action contribution S n is a generating function of the first kind as before, just not for a 
canonical transformation (as the corresponding Hamiltonian map is not symplectic but rather 
pre-symplectic), but for a singular global time evolution. 

Nevertheless, it will generally be convenient to work with suitable configuration and phase 
space extensions such that one obtains phase spaces of equal dimension at (n — 1) and n. This 
simplifies the discussion of the consequences for the preservation of the symplectic structure 
under H n -i- 

The prescription for the phase space extension is simple. Assume there exist 'old' configura- 
tion variables which appear at time step (n — 1) but do not have counterparts x° n at step 
n. We would like to extend the phase space at time n such as to include the pairs (x°,p™). To 
this end, extend the configuration manifold Q n to a new configuration space Q n := Q n x Q^f*, 
where Q^ 1 is a suitable configuration manifold of appropriate dimension coordinatized by the 
desired x° H In this way the coordinate form of the Lagrangian two-form (13. 5p will be again 

7 As will be explained later, the pre-constraints in V n and the post-constraints in V n +i each define a first 
class set of constraints. The pull back of the symplectic form to a constraint hypersurface described by first class 
constraints is a pre-symplectic form with one null direction per (independent) constraint. 

8 For instance, in Regge Calculus one has a priori, i.e., before imposing any equations of motion or generalized 
triangle inequalities, Q n — R+ x ■ ■ ■ x R+ (the lengths cannot be negative). One can, therefore, simply extend Q n 
by some additional products of R+ coordinatized by the lengths l„ of 'old edges' (which actually do not appear 
in the triangulation at step n). See [2] for more details. 



13 



a square matrix, however, with (additional) left and right null vectors corresponding to the 
formally added variables x n . The action S n as well as S n +i do not depend on these variables 
as these were only introduced for book keeping purposes. The x° n are thus both a priori and a 
posteriori free parameters A° = [x° n according to definition 13.31 Upon extension to Q n we are 
therefore in the case of singular discrete dynamics and all the results of subsections 13.3.11 and 
13.3.21 apply. The extension for 'new' variables x n £W which only appear at step n but do not have 
counterparts at (n — 1) proceeds analogously. 

On the extended configuration spaces Q, n -\ x Qn, where now dim Q n -i = dim Q n , one may 
then perform the discrete Legendre transformations (|3.10l 13. lip , yielding a mapping to the 
extended phase spaces V n -i '■= T*Q n ~i and V n := T*Q n . Denote by (a^f^Pnew) an< ^ i x n ,Po) 
the pairs by which the phase spaces at re - 1 and n have been extended, respectively. The 
Hamiltonian time evolution map 7i n -i is extended to the formally added variables by simply 
using the action S n as a generating function of the first kind with trivial dependence on x™^, x n , 
resulting in the evolution equations 

-n-1 dS n = „ + n _ dS n 

fnew Q_tifii]i u ' f' 



new <9x ne ™ fnew Q x new ' 

- n n-i _ dS n + n n -^IL-n (-12^ 

That is, additional pre- and post-constraints 

- r n-l _ - n n-l _ 9S n _ + _ + _ dS n _ 

'-'new - Pnew ~ ~ Q x n ' ~ U ' °o ~ Po ~ ~ U {(5.10) 

now arise at n — 1 and n, respectively, which are in number equal to the difference in dimension 
between the extended and unextended configuration spaces, respectively. The total resulting 
numbers of pre-constraints at n — 1 and post-constraints at n are identical as the Lagrangian 
two-form tt n (I3.5P corresponds to a square matrix which has equally many left and right null 
vectors L n ^i,R n (|3.22p . 

Thus, the added variables — being both apiori and a posteriori free — cannot be determined 
and their conjugate momenta are constrainedly We are free, however, to make the 'gauge choices' 
x neuj _ x new an( j x <^_ 1 = x °. Since S n does not depend on x^l^x", neither does any of the 
already existing pre- or post-constraints at n — 1 and n, respectively, such that both new sets 
in (|3.26p are first class. Hence, returning to the unextended ('reduced') phase spaces V n -\ and 
V n is simply performed by a partial reduction procedure consisting of imposing ~Pnew = an d 
+ Po = an d factoring out the corresponding first class flow. This amounts to simply dropping 
the added pairs (a^-liPneto) an d {% n iPo) from their respective phase spaces. 

In that sense, the phase space extension is a trivial extension which simply adds constrained 
canonical pairs. Nevertheless, precisely this allows us to directly apply the discussion of singular 
systems of the previous section to systems with varying number of variables. 

Example 3.2. Coming back to our massless scalar field, we found the time evolution equations 



9 If one also extended the phase spaces at previous time steps by the added variables x new , the momentum 
e-w would still be vanishing. The same holds for the variables x° and later time steps. 
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Figure 4: A global time evolution move for a 
scalar field on a triangulation with Q n = 2 and 
Qn+i = 3. Dashed lines represent periodic iden- 
tification. 



from step n to (n + 1) as 



7T 



7T; 



n+1 




ft) 



+ 



2 ( Pn+l 



b^ 1 



(3.27) 



with ^ij^ 1 encoding the adjacency relations between the vertices i at time n and the vertices 
j at time (n + 1). The number of configuration variables is Q n and Q n +i for times n, (n + 
1) respectively. Thus, if these numbers disagree we can introduce auxiliary additional field 
variables. Since the corresponding dynamics is trivial we will refrain from doing so explicitly. 
The Lagrangian two-form agrees with the adjacency matrix A?J~ . If Q n +i = Qn + N with 



N > we have at least N right null vectors (i? n +i)r 



1, . . . , N of the Lagrangian two-form 



(and N trivial left null vectors if N auxilliary variables <^>° are introduced). The corresponding 
iV post-constraints are given as 



+(jn+l 



E 



r n+l 




+ 1 



J n+1 



+ 



2^ n+1 



2Vn+l 



{Rn 



+1> 



(3.28) 



For instance, for the triangulation in figure H] with two vertices at time n and three vertices 
at time (n + 1) we have the adjacency matrix and right null vector 



.4 



n+1 



1 

2 1 1 



Rn-\ 




(3.29) 



This yields the post-constraint at time step (n + 1) 



C 



n+1 



TV, 



n+1 n+1 5 (±2 



7T 



(4>l +1 - 4> n+ i) , 



(3.30) 



corresponding to an a priori free parameter A n+ i = 4> 2 n+ i — 4> n +i- 

Time-reversing the triangulation, i.e. just exchanging n and (n + 1), one obtains a pre- 
constraint at time n 



- C n 



n „n i 5 / ±2 



7Tn — 7T 



+ |(0n-^), 



(3.31) 



corresponding to an a posteriori free parameter fi n - 



"n — Y n Yn- 
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S n _i Figure 5: A fat slicing of a triangulation. 



3.5 Bulk variables and Hamilton's principal function 

Some systems involve bulk variables, i.e. variables — which we shall denote by x l n — that only 
appear in a single action contribution S n . Thus, the corresponding equations of motion simply 
read = 0. Examples appear in the context of evolution schemes in Regge Calculus which, 
e.g., involve fat slices [2] or tent moves [21 (HJ HHJ [33] . For example, figure [5] depicts two adjacent 
fat slices bounded by non-intersecting hypersurfaces S n in a triangulation. The variables are 
given by the lengths of the edges and all the lengths of the dashed edges will all appear as bulk 
variables either in S n or in S n -\. 

Such bulk variables also appear when performing 'effective' evolution moves n — > n + k which 
consist of multiple intermediate basic time steps n + l,n + 2,...,n + k — 1. The variables at these 
intermediate time steps will only appear in the action S njn+ k = Yl?=n+i &l but not in the other 
action contributions, describing other (possibly effective) time evolution moves. Integrating out 
these variables would lead to Hamilton's principal function, which, in fact, is a function of the 
boundary variables only. In the following we shall relabel the time steps so that the effective 
evolution move appears as n — > (n + 1). 

By virtue of this splitting of the configuration variables into bulk variables x l n and boundary 
variables x^, which in general will appear in two action contributions, the evolution equations 
defining H n -\ (I3.13P become 



-vT 1 



d_Sn 





_ dS n 








dS n 


Pi 





~ vT := ~^t~ 1 =0. 

As a result, one obtains the pre-constraints ~C f n_1 = = (S n does not depend on 

which, by momentum matching, leads to + Pt ~ 1 = 0. 

Thus, also the equations of motion for the bulk variables are implemented by momentum 
matching, in this case by a constraint, 

o = - P r l = + P? = ^r- (3-33) 

This, in fact, is a general attribute of the present formalism: equations of motion manifest 
themselves as canonical constraints p n = 0. 

Denote by Xn( x n -ii x n > Km ) the solutions to the equations of motion (13.330 for the internal 
variables x l n where the K m denote possible free parameters. Substituting this into S n , one can 
define an effective action S n (or Hamilton's principal function) which only depends on the 'true 
boundary' variables x e n _ x and x e n . The effective action defined in this way does not depend on 
the free parameters n m : The solutions X n ( x n-ii x ni K?n ) extremize the action and the effective 
action is equal to the value of this extremum. The free parameters appear only if the extremum 
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is actually not a proper extremum, rather there are constant directions along which the action 
(extremum) does not change. Below, however, we will redefine the effective action for the case 
that the equations of motion for the bulk variables imply constraints on the boundary variables. 

The effective action defines the effective dynamics of the remaining degrees of freedom and, 
in analogy to definitions 13.11 and 13 . 21 the effective pre- and post-Legendre transformations ¥~S n : 
Sn-l x Qn — > T*Q n -i and ¥ + S n : Q n -i x Q„ -> T*Q n , respectively, where Q n denotes the 
configuration manifold at step n after solving (|3.33p . The following theorem shows that the 
dynamics defined by the effective action S n is equivalent to the dynamics defined by S n . 

Theorem 3.2. Let Q n := Q n x Q n , where Q n and Q n are coordinatized by x n and x n , respec- 
tively. The following two diagrams commute 

F+ F~ 
Sn-l x Q n > T*Q n Sn-l X Q n > T*Q n _! 



eom 13J3\ 



pj =_r pJ 



V v V \|/ 

v+ f- 

Sn-l x Q n > T*Q n Q n -! X Q n > T* Q n -l- 

This also holds in the presence of constraints on the boundary data, i.e. on Sn-l x Qn- 

Proof. To show that the effective action defines the same dynamics as the action we started with, 
we have to convince ourselves that the momenta and constraints agree if we perform either of 
the following two procedures: 

(a) We first define the canonical time evolution (|3.32p and impose there the momentum matching 
~Pt = = + Pt, which implements the equation of motions for the bulk variables x n (right-down 
in the diagrams). 

(b) We solve the equations of motion for the bulk variables and define the effective action. We 
use this effective action to define the canonical time evolution from (n — 1) ton (down-right in 
the diagrams). 

The equivalence of the two methods in the absence of constraints on Q n -i x Q n follows from 

-~n-l _ 8S n _ 0S n 0x n _ 9S n _ - n-l (o oa\ 

ux n-l UJj n-l ux n x n-l udj n-l 

where the last equation holds due to the equations of motion = for the bulk variables. A 
similar calculation shows equality of the post-momenta. 

However, this argument is not sufficient if the equations of motion for the bulk variables 
impose constraints between the boundary variables x n and x n _ v These can come as proper 
canonical constraints involving only configuration variables at one time — either x n or x n-1 . 
Such constraints are called holonomic and arise in this form after momentum matching for the 
bulk momenta. The other type of constraints that can occur are relations between variables at 
time (n — 1) and time n and are thus not proper canonical constraints. (Canonical constraints 
are equations of motions involving canonical data at only one time.) We will refer to both kinds 
of constraints as boundary data constraints. 

Let us assume that after solving for the bulk variables one finds M independent boundary 
data constraints H^x^^jX^) = (forming an irreducible set). The form of these constraints is 
not unique — but locally these constraints can be solved for M variables x h = x h { x n -i' x n) ou t 
of the x n _ 1 ,x n , where e denotes the remaining variables at (n — 1) and e the remaining variables 
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at n. We furthermore assume that the matrix of derivatives |^ is (locally) invertible. This 
allows us to express the derivative of x h ' with respect to x e n _ x 

dx h ' = fdHhy 1 dH h 
dx^_i \dx h ' J dx^_\ 

which can be derived by taking the derivative of Hf l (x^ l _ 1 , x^, X h { x n-n x n))- 

Also, if we have the equations of motion for the bulk variables x l n resulting into M (inde- 
pendent and irreducible) relations restricting the boundary data, we will have at least M free 
parameters K h (among the K m introduced below f|3.33|) ) for the solutions Xn- We will take these 
additional bulk variables for the effective action, acting as Lagrange multipliers. If 
there are more than M free parameters a suitable choice has to be made. Below we will derive 
a condition for this choice. 

Thus, for the effective action to impose the boundary data constraints we define 

Sn(x n _i, X n , K n ) := S n (x n _ l ,X n ,X (x n _i, X n ), Xn( x n-li x ni K )) + 

A- k ( K n 1 x n-li x n)Hh{ x n-li x n) ■ (3.36) 

The functions A h will be determined such that the momenta obtained by the two different 
methods (a) and (b) coincide. But we need to assume that Q^jr is invertible so that the equations 

of motion for the set {k^} impose the constraint set {H h }. (The first term on the right hand 
side of (|3.36p . as previously argued, does not depend on K m .) 

To see this, consider the momenta at (n — 1) as defined from the effective action (we omit 
terms proportional to Hh as these vanish if the equations of motion for the are imposed) 

ds n ds n d x h kh dH h 



6 d x n-l 3 xh d X n-l d X n-l 

dS n dH h ( kh dS n (dH h 



d x n-i ® x n-i \ d xh ' \dx h ' 
-pi}- 1 = -K h ^L- (3.37) 

where x\_ x are those variables among the x h which are associated to the time (n — 1). 
Comparing with the momenta obtained from the original action, 

- 9Sn ~„n-l _ 9S n ( ooo\ 

PS ~ dxU' H ~ dxti ( } 

(evaluated on the surface {x l n = Xn,H h = 0}), we see that we need to set 

to obtain for the pre-momenta p n ~ l = p 11 ^ 1 . The same line of arguments can be made for the 
post-momenta. 

Using the solutions x l n = Xn( x n-n x ni K?n ) an d xh = X h { x n-n x n) lli (|3.39|) will make the A h 
dependent on kJJ . The assumption on the invertibility of translates then into the assumption 
that 

d dS n \ dxn 



dx l n dx h ' } dn h 



(3.40) 
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is invertible. This, in particular, might restrict the choice of the M parameters K h among the 
free parameters K m for the solutions x l n = Xn- 

This finally shows that the dynamics defined by the original action and the effective action 
is equivalent. □ 

In the case of constraints on the boundary variables one has to be careful to incorporate 
these constraints into the effective action. A typical example for a bulk variable x l n is provided 
by Lagrange multipliers. A constraint on boundary data arises if we have an action of the form 

S n (x n _i, x n , x n ) = S n (x n _ 1 ,x n ) + x n H(x n _i,x n ). (3-41) 

In this case x l n will be undetermined by the equations of motion for x n . Thus, we will have 
x l n = k. For A we find according to ()3.39|) 

A " h+ dx^{dx^J ■ (3 ' 42) 

Hence, A replaces the bulk variable that acted as a Lagrange multiplier in the original action. 

Another example are effective evolution moves n — > n + k, where we treat all intermediate 
variables as bulk variables. The previous discussion implies that it does not matter 

(i) in which ordering one integrates out intermediate evolution steps, and 

(ii) in which ordering one engages Legendre transformations and equations of motion 
in order to evolve from n to n + k. 

Finally, in appendix^] we show that, under the simplifying assumption of absence of bound- 
ary data constraints, the effective canonical two-form w n on T*Q n and the effective Lagrangian 
two-forms (l n on Qn-i x Q n and f2 n +i on Q n x Qn+i are, indeed, correctly related by pull-back 

n n = (¥ + S n )*oo n , U n+1 = (¥~S n )*Co n (3.43) 

via the effective post- and pre-Legendre transformations F^Sn, respectively. 



4 Local dynamics of variational discrete systems 

In the previous section we have only considered global time evolution moves, i.e. time evolution 
moves such that no subsets of variables at any two neighbouring time steps n, n + 1 coincide 
(except possibly in the boundary). In this section, let us now also consider local evolution moves 
which only evolve subsets of the canonical data of a given time step and which for distinction 
from the global moves we shall label by k € Z, instead of n € Z. That is, different discrete 
time steps k, k + 1, may partially overlap such that they involve coinciding subsets of canonical 
variables. Such time evolution moves occur, e.g. in simplicial gravity or lattice field theory 
when only small regions of a discrete hypersurface E are evolved in discrete time such that 
Sfc n Xfc+i ^ 0; the Pachner moves [2] and the tent moves |18^ [2] are particular examples of 
such local moves (see these reference for their canonical formulation). Such local evolution 
moves define the most basic evolution steps from which more complicated time evolutions can 
be constructed. In particular, the Pachner moves constitute an elementary and ergodic class of 
local moves applicable to arbitrary triangulations such that one can map between any (finite) 
triangulations of fixed topology by finite sequences of these moves [Ml [35] . 
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4.1 Coinciding subsets of variables and momentum updating 

We will formulate the evolution equations corresponding to the local evolution moves directly 
in the Hamiltonian picture. The reason is twofold: 

(i) In the discrete the Lagrangian picture (see section [3. l.lj) necessitates the complete config- 
uration data of two consecutive time steps. But the action contribution S k+ i governing 
a local evolution move k — > k + 1 only contains data which is directly involved in the 
evolution move. 

(ii) Subsets of data in Q k and Qk+i coincide. Accordingly, the contribution S k+ i alone does 
not suffice to define the Legendre transformations (definitions 13.11 and 13. 2p from Q k x Q k+ \ 
to the phase spaces V k ,V k +i. 

In the Hamiltonian picture the configuration data of the second time step are replaced by the 
momenta which are defined at the same time step as the configuration data. That is, in the 
Hamiltonian picture we need only one time step to encode the canonical data; if subsets of 
variables coincide at consecutive steps, the canonical data merely need to be appropriately 
updated in the course of the evolution move (while the momenta themselves are defined using 
the global structure of section |3.1.2|) . 

In section 13 .1.21 we saw that the action contribution S n+ \ of a global evolution move n — > n+1 
is the generating function of the first kind of the corresponding global Hamiltonian time evolution 
map T~L n . In this section, we shall introduce the generating functions of the local Hamiltonian 
time evolution maps, which, for distinction, we shall denote by t) k : V k — > Vk+i, corresponding 
to local time evolution moves k — ¥ k + 1. 

Firstly, suppose some set of variables x b appears in several time steps, e.g., such that x k+l = 
x k . Assume x b is not involved in the local time evolution such that 

4 + i = 4, p k h +1 = pt (4.i) 

These evolution equations cannot be generated by the action, since S k+ i depends on neither 
x\ nor on x k+l , given that x b is not included in the local dynamics. Instead, the identity 
transformation — either governed by a generating function of the second (depending on old con- 
figuration and new momentum variables) or the third kind (depending on new configuration and 
old momentum variables) — is appropriate: 



nn(-r k r, k+1 \ - -r b n k+1 n k - — — - r, k+l <r k+1 - — - >r k 

< ~ r 2{X b ,p b ) — —X k p b , p b — — k — p b , X b — — k+l — X b 

ox b op b 



„ , k+i k\ k fc+i dGs k k dG?, b 

Pb = ~g^b = Pbi x b = -Q^p = x k- 



Next, let us consider the case in which some configuration variables do not evolve x k = x k+1 , 
yet in which either x k or x k+1 appear in Sk+i (but not both for the same index e) such that 
their conjugate momenta are transformed. The additivity of the action implies either 

n k - n k+i _ 9S k+1 (x k ) k+1 _ k dS k+1 (x k+1 ) 

Pe ~ Pe g x e ™ Pe ~ Pe + > W 



which we shall call (|4.2j) momentum updating. The last definitions (|4.2p can be summarized into 

Jfc+i =n k , dS k+1 {x k+1 ) dS k+l {x k ) 

Pe Pe ' o e a e " ^•°' ) 

ux k+l ux k 
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One of the derivatives will be zero, the remaining terms define both ways of momentum updating. 

The generating function for momentum updating is of second or third kind: either add Gi 
or G3, respectively, to Sk+i, the latter of which either depends only on the old configuration 
variables x| or only on the new configuration variables The same construction also applies 

to the case of varying numbers of canonical variables when formulated on extended phase spaces, 
as we shall see in the following subsection. 

4.2 Preservation of symplectic structures and constraints under momentum 
updating 

The preservation of the symplectic structure, as proven in theorem 13. II for singular systems and 
discussed in section [3^1 for varying phase spaces, holds for global evolution moves where variables 
at n and (n + 1) do not coincide. We still have to investigate the preservation of the symplectic 
structure under local evolution moves which involve coinciding subsets of variables and proceed 
by momentum updating for all canonical pairs on the extended phase space. To circumvent the 
problem of changing phase space dimensions, we shall work on extended phase spaces Vkj'Pk+i 
as introduced in section 13.41 

At this stage we shall distinguish four types of local evolution moves. Any other conceivable 
kind of local moves can be treated in complete analogy!^ Let us detail the corresponding 
evolution equations case by case and subsequently consider the symplectic structure. 



Type I: The move introduces 'new variables' but does not remove 'old variables.' Assume 
that K 'new variables' arise. Accordingly, extend the phase space at time k by K pairs {xi,p^) 
which correspond to the K canonical pairs at time k+1. Additionally, pairs (x b ,pb) 

occur which do not change during this evolution move and variables (x e ,p e ) for which only the 
momenta are updated. The local Hamiltonian evolution map fjfc for type I is thus given by the 
following momentum updating 

4 = 4+1 , P$ +1 = p\ , (4.4) 

r e - r e w fc+1 - n k + dSk +^ X k+V x k+l) u ,s 

x k — x k+l > Pe — Pe ^ JT^ J l 4 -°J 

UX k+l 

Pi = 0, pT = ^±^kll^A . (4.6) 

ox k+l 

We choose Sk+i to be a function of the variables at time k + 1. Equations (|4.6p contain K 

constraints = p^ which are simultaneously pre- and post-constraints, as well as K post- 

dSi (x e x n ) 

constraints + C^ +1 = p^ +1 ' +1 g x n +1 ' k+1 ■ Both x^,x^ +1 are undetermined. 

Example 4.1. Consider a scalar field living on the vertices of a 2D space-time triangulation. 
The evolution move given by the 1-2 Pachner move (see figures [6l [7]) introduces one new vertex 
v and thus one 'new' field variable at k+1, while preserving all other vertices of Hence, 
this move is of type I. Using (|2.4p . the momentum updating map reads 

4 = 4+1, = 4, 6 = 1,4,5, 

4>t = 01+1 , ^e +1 = < + <Pt + i ~ \ {K+i + O , e = 2, 3 , 



7T„, = U , 7L, 



r-k _ n ^k+1 _ a,v L (j2 1 ^3 

10 



4 + 4>k+i) 



In particular, we shall ignore the case where 'bulk' variables are involved in the local move. These cases can 
be treated via the recipe provided in section 13.51 
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where e + l = 3ife = 2 and e + l = 2ife = 3 (see figure E] for further notation). 

Furthermore, the 1-3 Pachner move in 3D Regge Calculus and the 1-4 and 2-3 Pachner 
moves in 4D Regge Calculus are of type I (see [2] for details). 





fc+i 



■>k+i 



(a) 



(b) 



Figure 6: The 1-2 Pachner move glues a triangle onto the ID hypersurface Efc and introduces one new 
field variable </) k+1 at vertex v at step k + 1 and is of type I (see figure [7] for the ID perspective). 

Type II: The move removes 'old variables' but does not introduce 'new variables'. This is 
the time reverse of type I. If K old variables are removed, extend the phase space at step k + 1 
by K pairs (x k+l ,p^ +l ). We choose S^+i to be a function of the variables of time k. The local 
Hamiltonian evolution or momentum updating map f)fc for type II reads 



,J> 
•' k 



• r A- 



P k a +1 



x k+l > 



b k+l ) 



Pc 



Pb 



p: 



k+l 



dS k+1 (x e k ,x° k ) 



dS k+1 (x 



dx e k 

ki X k) 



dx° k 



(4.7) 
(4.8) 

(4.9) 



Equations (|4.9p contain K pre-constraints = p k H k+ Q^, k ' Xk and X constraints + Cq +1 = 

Po +1 > w hich are simultaneously pre- and post-constraints, while x k+l remains undetermined. 
We shall denote the partial pre-constraint surface defined only by the ~ C„ in the extended phase 
space Vk by KT C Vk- No further pre-constraints are created in the course of the move I 11 ! 

Example 4.2. Consider, again, the scalar field living on the vertices of a 2D triangulation. The 
2-1 Pachner evolution move (see figures [8] and [7|) is the time reverse of the 1-2 Pachner move 
and thus of type II. It removes a vertex v* and, accordingly, an 'old' field variable (jf k from S. 
Employing (|2.4p yields its momentum updating map 



7T 



4 



k+l 



^k+l i 



k+l 



Kb 



b 



0. 



7T, 



TT. 



k+l 



+ 



1,4 
2,3 



1 

+ 2 



Subsequent evolution moves may lead to additional pre-constraints on the same data. 



Figure 7: The 1-2 and 2-1 Pachner 
moves in ID. 
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where e + l = 3ife = 2 and e + l = 2ife = 3 (see figure [5] for further notation). 

Moreover, the 3-1 Pachner move in 3D Regge Calculus and the 3-2 and 4-1 Pachner moves 
in 4D Regge Calculus are of type II (see [2] for details). 





k+l 



(b) 



Figure 8: The 2-1 Pachner move glues a triangle onto the ID hypcrsurface such that vertex v* 
disappears at step fc+1 together with the variable <j>\ and is of type II (sec figurc[7]for the ID perspective). 

Type III: The move removes K 'old variables' and also introduces K 'new variables'. Ex- 
tend the initial phase space by K pairs {x k ,p k ) and the final phase space by K pairs (x k+1 ,p k+l ). 
The momentum updating map \] k is given by 



J> 

a k 


— x k+l ) 


P k b +1 


x\ 


= x k+l > 


P k e +1 


Pn 


= o, 


P L n +1 


,k+l 




= o, 


Po 



Pb 



Pe + 



dS k+1 (x 



fc+i' x fc' x fc-i-i/ 



dx\ 



J k+l 

dSk+\{x e k+v x° k , x k+l ) 

dx t+i 
dS k+ i (x e k+1 , x° k , x k+1 



dx° k 



(4.10) 
(4.11) 

(4.12) 

(4.13) 



(We could equally well have chosen to let S k+ \ depend on x k , instead of x e k+1 .) The p k n = 
and the p^ +l = each constitute K constraints that result from phase space extensions and are 

each simultaneously pre- and post-constraints. Furthermore, if the K x K matrix a oJ't 1 is of 

ax k ox k+i 

rank K — k, then k non-trivial pre-constraints at k and k non-trivial post-constraints at k + 1 
will arise from the second equations in (|4.12p and (|4.13p . respectively. Among the x k+1 , K — k 
variables will be determined via the second equation in (|4.13j) . while the remaining k variables 
will remain (a priori) free. Furthermore, the x k and x k+1 remain undetermined and can be 

The above k non-trivial 



x k+l an d x °k+l 
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gauge fixed to arbitrary values, for instance x k 
pre-constraints at k alone define the partial pre-constraint surface fC k C V k . 

Example 4.3. Consider a scalar field living on a 2D space-time quadrangulation where each Y> k 
is given by a 'zig-zag line'. Adding a square as a fundamental building block annihilates one 'old' 
vertex v* , introduces a 'new' vertex v and preserves the remaining vertices in the move k — > k + l 
(see figure [9]). Hence, this evolution move annihilates one field variable introduces one 

new field variable <f> k and is of type III. The action of a scalar field on a (rectangular) square 
reads 
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where i labels the vertices of the square and i + 1 = 5 coincides with i = 1 due to periodicity. 
The corresponding momentum updating map is thus given by 



4 






= > 


6 = 1,4,5 
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= x k+l > 


7T 

e 


= 7r* + 20g +1 -$f-$ +1 , 


e = 2,3, 


_k 
ir v 


= o, 


7T 


= - </>fc + l — 4>\ + l ! 




.fc+i 

V* 


= o, 


J* 


= -2$f + ^ + ^ 





(see figure [9] for further notation). Notice that + C k+1 := ir k+l — 2^ +1 + <ji>| +1 + 4>\ +l defines a 
post-constraint at k + 1, while "C** := ir k * + 2<^j£ — 0| — 4>\ constitutes a pre~constraint at k. 

Furthermore, the 2-2 Pachner move in 3D Regge Calculus and the tent moves (after solving 
the tent pole equation) in any dimension are of type III (see [21 [18] for details) . 




(a) (b) 

Figure 9: Gluing a square onto the ID hypcrsurface Efc annihilates a vertex v* and its field variable §" k 
and introduces one new vertex v and its field variable (j^k+i- This move is of type III. 



Type IV: No 'new variable' is introduced and no 'old variable' is removed in the course of 
the move. The momentum updating map t) k simply reads 

4 = x b k+1 , P k b +1 =P k b , (4.14) 

r e - r e r> k+l - n k + dSk+1 ^ X k+^ (a i z\ 

x k — x k+l > Pe — Pe H ■ l 4 - i;3 J 

(Sk+i may equally well depend on x k , rather than x e k+1 .) No pre- or post-constraints arise. 

Example 4.4. The 2-2 Pachner evolution move for a scalar field living on the vertices of a 3D 
triangulation is of type IV. No vertex is annihilated or introduced so neither is any field variable 
(see figure fTOj) . We abstain from detailing the corresponding momentum updating map. 



The following theorem shows that momentum updating preserves the post-constraints and 
the symplectic structure restricted to the post-constraint surface in the move k — ¥ k + 1. How- 
ever, for types II and III this constraint surface at k is crucially further reduced by intersec- 
tion with the partial pre-constraint surface /Cj7 El Hence, momentum updating — as presently 

12 For types I and IV PI = because no pre-constraints independent of the post-constraints arise in 
the move. In particular, for type IV = 7\. 
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(a) (b) 

Figure 10: Performing a 2-2 Pachner move in the 2D hypersurface E corresponds to gluing a tetrahedron 
onto the hypersurface as shown. This move neither removes nor introduces any vertices or field variables 
and is of type IV. However, the momenta conjugate to e = 1, . . . , 4, must be udpated. 



formulated — is a pre-symplectic transformation. In fact, there also exists an alternative way of 
formulating momentum updating as a canonical transformation on the full extended phase space 
which we shall briefly introduce in appendix [Bl 

Theorem 4.1. Let uok^k+i be the symplectic forms on the (extended) phase spaces Pk,Pk+i and 
let C^,C^ +1 be the post-constraint surfaces at steps k and k + 1, respectively. The momentum 
updating map h/% preserves the symplectic structure as follows 

where for 

Types I and IV: : — > C^ +1 and tk/k+i '■ ^t/k+i ^ ^k/k+i are embedding maps. 

Types II and III: : n K.^ — > C^ +1 and JC^ C Pk is the partial pre-constraint surface at 

k and tk '■ D K~j^ Pk an d tk+i '■ C^+i ^ ^k+i ar z embedding maps. 

In particular, maps all post-constraints on Pk to post-constraints on Pk+i- 

Proof. The proof is given in appendix O □ 

Consequence: The preceding theorem has important repercussions for the preservation of 
the rank of the symplectic structure and the number of constraints throughout the entire discrete 
evolution. While the local evolution moves of type I and IV preserve the symplectic structure 
restricted to the post-constraint surfaces before and after the move, local moves of type II and 
III only preserve the symplectic structure further reduced by the partial pre-constraint surface 
K,Z ■ This, in particular, means that type II and III moves will, in general, not preserve the post- 
constraint surfaces dimension. Rather, the number of post-constraints at k + 1 can only be equal 
or higher than the number of post-constraints at fcO the number is equal if the pre-constraints 
defining KT are such that they coincide with post-constraints at k or if they are of second class 
together with the post-constraints at k (second class constraints do not further reduce the rank 
of the symplectic form p3]). The number increases if some of the pre-constraints defining /C7 
do not coincide with post-constraints at k and neither are rendered second class. As a result, on 

13 The 2-2 Pachner move in 3D Regge Calculus — which is of type III — actually does not lead to any non-trivial 
pre-constraints at k (see [2] for details). That is, k — for this move and ZC^T = Vk such that the number of 
post-constraints at k and k + 1 are identical. 
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an evolving slice, the number of post-constraints will either remain constant or grow, but cannot 
decrease. 

This also has consequences for the number of constraints at some fixed global step n. Theo- 
rem [3J] showed that a global evolution between n — 1 and n preserves the symplectic structures 
restricted to C~_ 1 and because the number of pre- and post-constraints coincided. In con- 
junction with the present theorem 14. 1\ it implies that, if step n is evolved forward by a local 
evolution move k — > k + 1, where the local step k now coincides with the global step n, the sym- 
plectic structure restricted to the pre-constraint surface at (n — 1) and the new post-constraint 
surface at k + 1 must again be preserved. However, if the local move was of type II or III, it 
may have increased the number of post-constraints at k + 1 as compared to n = k. Conse- 
quently, also the number of pre-constraints at (n — 1) must actually have increased through the 
evolution from n = k to k + 1. That is, in general the number of constraints at fixed k can 
only stay constant or increase by further evolution. This is a consequence of solving equations 
of motion in the course of the moves; solving more equations of motion can lead to more, but 
not less constraints because the equations of motion can act as secondary constraints. This has 
important ramifications for the notion of propagating degrees of freedom and the reduced phase 
space. We shall discuss all of this amply in the next section [5j 



5 Constraint analysis in the discrete 



A discrete time evolution proceeds in discrete steps and therefore, in contrast to the continuum, 
is not generated by a set of constraints via a Poisson bracket structure which necessarily has an 
infinitesimal action. Rather, it is the time evolution moves which generate the discrete evolution. 
The exception are theories in which the diffeomorphism symmetry of the continuum is preserved 
as a vertex translation symmetry [7J EJ El EH] (and one considers a time evolution which preserves 
the connectivity of the lattice). This is the case for 3D Regge gravity and also for 4D Rcgge 
gravity at vertices embedded in a flat geometry. In general, however, the symmetry is broken 
in 4D Regge gravity [8] such that continuous time evolution has to be replaced by discrete 
steps, as also proposed in the consistent discretization program [231 [Ml [36]. Understanding the 
complicated status of constraints and symmetries in discrete 4D Regge gravity is one of the 
motivations for the present work. 

Apart from the issue of generating time evolution, the constraints in the discrete should 
assume similar roles to those of their continuum analogues [13], namely: 



role of constraints in the continuum 

(i) guarantee correct dynamics 

(ii) generate symmetries 

(iii) classify degrees of freedom 

(iv) generate time evolution 



role of constraints in the discrete 

(i) guarantee correct dynamics 

(ii) generate symmetries 

(iii) classify degrees of freedom 



It is the goal of the present section to show that this is, indeed, the case. To this end, a 
constraint analysis for variational discrete systems, analogous to the continuum Dirac procedure 
[HI E3], needs to be developed. We shall discuss roles (i)-(iii) of the constraints in the discrete 
below after first considering the preservation of constraints by the discrete evolution. For this 
purpose it is necessary to discuss global evolution moves, but towards the end we shall also refer 
to local ones. Since we describe systems with evolving phase spaces on extended phase spaces of 
equal dimension, no generality is lost by restricting to singular systems where dim Q n = Q Vn. 
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5.1 Preservation of the constraints 



In the previous sections we have seen that a priori free parameters A n (or a posteriori free 
parameters ix n ) and, consequently, arbitrariness in the canonical evolution arise in the presence of 
post-constraints (or pre-constraints) . However, some of this a priori (a posteriori) arbitrariness 
may get fixed a posteriori (a priori) by the condition of preservation of the constraints. 

Continuum 

Firstly recall that in the continuum the free Lagrange multipliers X m of some primary con- 
straints 4> m may become fixed by the condition of preservation of the constraints under evolution, 
which reads 4> m = {0 m ,-ff + X m (f> m >} — 0, where H is the (non-vanishing) true Hamiltonian 
of the system and ~ refers to the fact that this equation needs only to hold weakly, i.e. on the 
constraint hypersurface |13| . This condition can 

(a) be automatically satisfied, in which case no new condition arises and the \ m remain free, 

(b) lead to secondary constraints which are independent of the 4> m and A m , or 

(c) lead to restrictions on the X m in which case some of the a priori free parameters get fixed. 

Case (c) is only possible in the presence of second class constraints. A similar, yet slightly differ- 
ent situation arises in the discrete where time evolution is not generated by a total Hamiltonian 
and one has to cope with two constraint surfaces in a given phase space. 

Discrete 

Since time evolution in the discrete is generated by the time evolution moves and not via 
the Poisson structure, we need to discuss the behaviour of the constraints under the former in 
order to discuss their preservation. 

For translation invariant system^] one requires preservation of the constraints under time 
evolution. In this case the primary pre/post-constraints are the same for every time step, but 
it is possible that C~ 7^ C+. One has to make sure that the image of the total constraint 
hypersurface, namely the intersection of the post- with the pre-constraint hypersurface at time 
n, C n := C~ HC+, is a subset of the (same) constraint hypersurface at time (n + 1), i.e. that the 
constraints are preserved under time evolution. Otherwise, one has to add further (secondary) 
constraints, which for translation invariant systems again will be the same at each time step. A 
related discussion of constraints in these special systems has been given in [24] in the context of 
consistent discretizations. 

For a non-translation invariant system, a preservation of constraints in the above narrow 
sense does not need to hold. Consider a (bare) evolution move (n — 1) — > n and another (bare) 
evolution move n — > (n + 1). Recall that Ti n -i '■ C~_ 1 — > and 7i n : C~ — > C^ +l . When 
considering the composition of the two moves, (n — 1) — > n — > (n + 1), we have to solve the 
equations of motion, i.e. match the momenta, at n. Consequently, at step n the pre-constraints 
now also have to hold for the post-momenta and the post-constraints have to hold for the pre- 
momenta. That is, the image of the pre-constraint hypersurface C~ 1 under time evolution "H n -i 
must be a subset of the total constraint hypersurface C n := C~ PlC^ at time step n. If this is not 
the case, secondary pre-constraints have to be added at step (n — 1) to appropriately restrict the 
data at (n — 1) to map under H n -\ only into C n because not all canonical solutions (n — 1) — > n 

14 By translation invariant systems we mean systems governed by an action such that S„(x n -i, x n ) is in form 
identical for all n. 
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can be extended to solutions for (re — 1) — > n — > (n + 1)- These secondary pre-constraints at 
time step (n — 1) arise from the pre-image under time evolution 7i n -i of the pre-constraint 
surface C~ Likewise, for consistency of the evolution (n — 1) — > n — > (n + 1), the pre-image of 
the post-constraint surface C^ +1 under T-L n must be a subset of the total constraint hypersurface 
C n at n. Otherwise, secondary post-constraints at the final step (n + 1) arise from the image 
under time evolution 7~L n of the post-constraint hypersurface C+. Consequently, constraints can 
'propagate' forward and backward in discrete time (see also the examples in section |5T2|) . 

For a larger sequence of evolution moves (for which in general both dimC~ 7^ dimC+, as well 
as dimC n 7^ dimC n +i are possible), this process is repeated until all constraint hypersurfaces are 
stable or an inconsistency has been reached. For translation invariant systems, this prescription 
is equivalent to the preservation of constraints given in [24] , However, in contrast to the trans- 
lationally invariant case, the secondary constraints are, in general, not the same at each time 
step and, moreover, crucially depend on the choice of initial n' = i and final n" = f time step. 
For instance, onsider an effective evolution move i — >■ n described by Si n and another effective 
evolution move n — > / described by S n f, where Si n ,S n f are effective actions. All constraints 
at n — including possibly propagated constraints, that would arise via the elementary evolution 
moves contained in both i — > n and n — > f — will arise as the post-constraints of Si n and as 
the pre-constraints of S n f. That is, the effective actions automatically contain all propagated 
constraints. 

Therefore, if we are using effective actions, consistency of the evolution simply requires to 
take the intersection of the post-constrained surface defined by Si n with the pre-constraint 
hypersurface defined by S n f as the total constraint hypersurface at time n, C n := C~ PI C+, 
for any elementary or effective evolution i — > n — > f (see figure [TO) . In addition, secondary 
pre-constraints at the initial step i may arise as the pre-image under time evolution of C n 
and, likewise, secondary post-constraints at the final step / may arise as the image under time 
evolution of C n . If, however, we now were to evolve further to i' < i and/or /' > /, additional 
propagated secondary constraints at any of the time steps may appear. 

As regards the canonical data at n, the imposition of pre-constraints in addition to the 
post-constraints leads to conditions which either 

(a) are automatically satisfied (i.e. the pre-constraints are dependent on the post-constraints), 

(b) are not automatically satisfied, yet which do not fix the flows and a priori free parameters 
A n of the post-constraints, 

(c) fix some a priori free parameters X n (likewise for the a posteriori free parameters /i n ), or 

(d) cannot be simultaneously satisfied such that C n = and the dynamics is inconsistent. 

This enforces restrictions on the discrete time evolution and the amount of arbitrariness in the 
data. The consequences of cases (a)-(c) for the dynamics, symmetries and observables, as well 
as roles (i)-(iii) of the constraints shall be discussed in detail in the remainder of this section. 

5.2 Restricting the dynamics 

Let us begin by briefly discussing role (i) of the constraints: ensuring the 'correct dynamics'. If 
we ensure that both the pre— and post-constraints (including possibly propagated constraints) 
are satisfied at each step, we obtain the correct dynamics because by momentum matching all 
equations of motion will be implemented. 

They arise from the pre-image of C~ because the pre-image under H n -i of C„ is by construction already the 
primary pre-constraint surface C~_ 1 . 
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(a) 



(b) 



Figure 11: (a) The •pre- constraint surface C~ and the post-constraint surface in the phase space 
T*Q n , in general, do not coincide, i.e. C~ ^ C+. In order to ensure the correct dynamics, we have to 
impose both the pre- and post-constraints at step n and thus must restrict to the intersection C n = 
C~ n C+. (If C n = 0, the dynamics is inconsistent.) (b) Matching symplectic structures at step n. 

In non-translation invariant systems the numbers of pre- or post-constraints at fixed n, in 
general, depends on the initial and final step between which one evolves: secondary constraints 
at n may arise as a consequence of imposing equations of motion at neighbouring steps. This, 
in fact, can only happen if case (b) above occurs at neighbouring steps. 

Example 5.1. (A) Consider the massless scalar field evolving in two moves as depicted in 
figure [121 (a). There are three scalar field variables at time n = 0, three at time n = 1 and two 
scalar field variables at the final time n = 2. The second move from n = 1 to n = 2 corresponds 
to the time reversed version we discussed in example 13.21 in section 13.41 Thus, we have one 
pre-constraint at time n = 1 



which leads to the free parameter = <f>\ — cj)\. For the evolution from time n = to n = 1 we 
have to consider the adjacency matrix 



which is invertible. Therefore, no pre- or post-constraints arise from the evolution — s> 1. The 
pre-constraint at time n = 1 has to be propagated back to a (pre-) constraint at time n = 0, 
which will also result into a free parameter fiQ, whose value cannot be postdicted by the canonical 
data at time n = 2. The so obtained constraint at n = will coincide with the pre-constraint 
of the effective time evolution from n = to n = 2, as illustrated in figure [T3l 

(B) On the other hand, consider the example in figure [12] (b). Here we replace the time 
evolution move — > 1 from the previous example by the time reversed time evolution move 
1—^2. From the considerations of the single moves in example 13.21 of section 13.41 one obtains 
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both a pre-constraint and a post-constraint at time n = 1 

-C 1 = wl-nl 

+ C i = ^ -4-1 - 4) . (5.3) 

These constraints are second class with respect to each other and indeed fix the a priori and a 
posteriori parameter Ai = fit = 4>\ — 4>\ = 0. The constraint hypersurface can equivalently be 
described by the two constraints C\ = 4>\ — 4>\ an d C\ = tt^ — tt^. Propagating the constraints 
forward to n = 2 and backward to n = does not lead to any constraints at these time steps. 
The effective time evolution —> 2 is thus regular. 

(C) Having the same number of variables at final and initial time steps does not necessarily 
guarantee a regular time evolution. If we exchange the two time evolution moves of the previous 
example (see figure [T2l (c)). we will obtain a pre-constraint at time n = and a post-constraint 
at time n = 2 

-c° = *§ - *§ + § (*§ - *g) , 

+ C 3 = vrf - vrf - | (<% - 4) (5.4) 
with the corresponding free parameters and A3. 




(a) (b) (c) 

Figure 12: Schematic illustration of (a) example 15.11 (A), (b) example 15.11 (B). and (c) example l5.ll (C). 
The numbers label the vertices and dashed lines indicate periodic identification. 

One might be surprised about the feature that constraints propagate backwards in time 
evolution. Therefore, let us provide an intuitive picture: we can imagine our evolution on the 
cylinder also as a radial evolution on an annulus region in two-dimensional space. If we proceed 
from time slices with more vertices to time slices with less vertices, we can interpret this as 
evolving from the outer boundary to the inner boundary of the annulus. Proceeding further 
and further with the evolution, more and more equations of motions in the ball defined by the 
outer boundary of the annulus have to be satisfied. Thus, we have to expect more and more 
constraints fixing the momenta as function of the boundary scalar fields. Closing the inner 
boundary we expect a totally constrained system at the outer boundary because the values of 
the scalar field on the spherical boundary should be sufficient to specify a solution uniquely. 
On the other hand, all the fields at the outer boundary correspond to free parameters since 
the fields at the outer boundary cannot be determined from the zero-dimensional phase space 
which arises after one has closed the annulus to a sphere. 
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(a) The two individual 
evolution moves — > 1 
and 1—^2. 



(b) Matching the symplec- 
tic structures at n — 1. 



(c) An effective evolution 
move — > 2 results with 
a new pre-constraint sur- 
face CrT at n = 0. 



Figure 13: Schematic illustration of cxample l5.ll (A). 



These examples highlight how constraints can severely restrict the space of solutions and 
how additional constraints at fixed n can arise. 

Now consider a general evolution up to step n. Assume a new evolution move n — > n + 1 is 
to be performed and one finds that a pre-constraint arising in this move is in conflict with the 
underlying canonical data at n. Then there are four options: 

1. Accept that we cannot perform this evolution move. In this case, perform some other 
evolution move 1^1 

2. (a) Change the underlying data by varying parameters which are a priori free up to step 

n such that the attempted move is possible, otherwise 
(b) restrict the space of initial data leading to n such that the attempted move becomes 
possible. 

3. Neither of 1 or 2 is possible and the evolution becomes inconsistent and stops. 

Option 1 is what one would choose if one solved an initial value problem, while options 2 (a) and 
(b) are what one would choose in case one attempted to solve some boundary value problem. 

5.3 Constraints and symmetries 

Next, let us elaborate on role (ii) of the constraints in the discrete, namely, generating gauge 
symmetries (if present). To begin with, consider the Hessian of the action, which is the matrix 
of second derivatives of the action with respect to 'bulk' variables. At time step n it reads 

= JPSn_ + ^H±L . (5.5) 

dx 3 n dx l n dxndx % n 

Null vectors of the Hessian define flat directions at the extrema of the action and thereby define 
genuine gauge directions. The Hessian thus plays a central role in the determination of gauge 
symmetry. Here we consider the Hessian with boundary data at (n — 1) and (n + 1) fixed. One 
should be aware that under further evolution, i.e. treating these boundary data as dynamical, 
the Hessian becomes a larger matrix (including blocks corresponding to the 'bare' Lagrangian 



16 For instance, a type I local move (in particular, in Regge Calculus a 1-D Pachner move, where D > 2 is the 
dimension of the spacetime triangulation [2]) is always possible because no pre-constraints arise. 
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two-forms at n and (n + 1), etc.), and the previous null vectors may cease to be null vectors of 
the effective Hessian. 

For non-linear theories, null vectors may only exist for the Hessian evaluated on solutions 
and the space of solutions may change if further evolution is taken into account (as a result of 
propagated constraints). Even more so, null vectors may exist only for special solutions. This 
is actually the case in 4D Regge Calculus, where null vectors exist for flat but not for curved 
solutions [8]. In the specification of irregular systems (|3.22p . however, we assumed a constant 
rank of the Lagrangian two-form in an open neighbourhood in the space of initial data. We 
shall, therefore, continue to make the assumption that the rank of the Hessians and Lagrangian 
two-forms is locally constant. Nevertheless, all considerations remain valid, even if a situation 
as in 4D Regge Calculus occurs, if one linearizes the theory around a solution (with additional 
symmetries) . 

Recall that an a priori free parameter A n and an a posteriori free parameter \i n were defined 
as functions of the configuration data at time step n, whose values cannot be pre- or postdicted 
from the canonical data at (n — 1) or (n + 1), respectively. A gauge models at time n is given 
by a free parameter X n = fj, n which is both a priori and a posteriori free and that can never be 
pre- or postdicted by any other data. 

Thus, given a solution s satisfying the equations of motion 

9(S n -i + S n ) d(S n + S n+ i) d(S n+ i + SV1+2) „ „x 

a = °> a = °> a = °' ( 5 - 6 ) 

dx n _i dx n ox n+ i 

the free parameter describes an infinitesimally displaced solution (s) n + eV n , V n £ TQ n at time 
n, which also satisfies the equations of motion. Taking the derivatives of the equations of motion 
in the direction of V n , we see that to each gauge mode there corresponds a vector V n at step n, 
which is a (i) null vector of the Hessian at n, (ii) a right null vector of the Lagrangian two-form 
f2 n , and (iii) a left null vector of the Lagrangian two-form 0, n +\. 

Next, let us study the relation between the presence of gauge symmetry and the first and 
second class nature of the pre- and post-constraints. To this end, we label the right null 
vectors {R n ) l r of £l n and their corresponding post-constraints by r = 1, . . . , and the iV£ left 
null vectors (L n )l of fJ n +i and their corresponding pre-constraints by I = 1, . . . , iV£. 

Theorem 5.1. The set of pre-constraints at fixed n and the set of post-constraints at fixed n 
each form a first class Poisson sub-algebra 

{-c?,-cp} = 0, {+c r va™} = 0, vz,/vy, (5.7) 

where the first equations holds on the pre-constraint surface C~ and the second on the post- 
constraint surface . 
Furthermore, 

{"CP, + d ~ 7l l '(L n )i, H% p/(R n )i, , W,r , 

where ji l (x n ,p n ) and p r r (x n ,p n ) are coefficient functions determined by the gradients ofC~ and 
C+, respectively and ~ denotes a weak equation valid on the constraint hypersurface C~ 



17 generally only denned in the linearized theory 
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Proof. The iV£ pre-constraints CV 1 and the post-constraints + C™ are defined through the 
equations 



o = ~c?(x n rp n ) 



-—5——{Xri ) XT l +l) 



= +C?(x n +p n ) (5.8) 

+ P"=§f^n-i,Xn) 

for arbitrary x n _i,x n and x n+ i. Differentiating the first equations with respect to x n+ \ and x n , 
respectively, and the second one with respect to x n _i and x n , respectively, yields 

o = arc? d 2 s n+1 _ arc? d-cf d 2 s n+1 
d~p] dxidx^ ' dxi d-pf dxidxi ' 

d+c? d 2 s n n 9+c r n , d+C? d 2 S n 



Note that the right hand equations also imply non-vanishing gradients g- p L and g+ p n 



because otherwise the constraints would be constants. These are equations between functions on 
Qn x Qn+i and Q n -i x Q n , respectively. We used the pre- and post-Legendre transformations 
in (15.8|) to map these configuration spaces to the phase space V n . The following equations, 
therefore, hold on the image of these maps, i.e. on the pre- and post-constraints surfaces. 

The equations on the left hand side of (|5.9l 15. lOf) specify null vectors of the Lagrangian 
two-forms, for which we earlier assumed a basis. Hence, 

r)~C n r)~^r n 

^-L = 7 /'(x n ,p")(L n )f„ — j, = p/(x n , P n(Rn%, (5.11) 
u Pj c„ u Pj c„ 

where 7/', p r r ' are appropriately chosen matrix coefficient functions. This associates a (non- 
vanishing) left null vector of f2 n +i to every pre-constraint and a (non-vanishing) right null 
vector of Q n to every post constraint. 

Consider the pre-constraints at n. Using (|5.9l I5.ll]) . we directly compute 



dx l n dp? dp? dx % n 



C n OXnOX l n OX\ l OX J n 

= 0. 

Likewise, using (15.101 I5TTT1) . one finds that the Poisson brackets between post-constraints at n 
vanish on the post-constraint hypersurface C+. 
Analogously, by (J51fll5.il]) . 

f)~r ,n rt+r 171 f>~r' n f)+c in 
s-nn +nn\ o L, l o u r _ o u ; o u r 

1 ' ' r } dxi dp? dp? dx\ 

_ l ' (r y I d 2 Sn+l _|_ 9 2 S n \ r' \j /g -j^) 

c^nc+ n 1 ydxidx^ dxidx^J 

where the term in brackets defines the Hessian of the action. This proves the statement. □ 
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In consequence, 

(i) a pre-constraint ~C™ is first class (in particular, its Poisson brackets with all post- 

ed" C n 

constraints vanish on the constraint hypersurface) if q- L ~ 0, 

(ii) a post-constraint + C™ Poisson is first class if Hfj d+p n+\ — 0. 

It is also possible that a pre-constraint is first class despite the Hessian not having null 
vectors. This happens if the corresponding left null vector 7/ (L n )\, is orthogonal to all right 
null vectors (R n ) l r i with respect to the metric defined by the Hessian. 

There is an immediate consequence of theorem 15.11 (applying to case (a) of section |5~T]) : 

Corollary 5.1. Let ~C n be a pre-constraint and + C n be a post-constraint. If this pre- and 
post- constraint coincide, i.e. C n := ~C n = + C n , then C n is necessarily first class. 

The following theorem implies that such coinciding pre- and post-constraints generate gauge 
symmetries of the action (or Hamilton's principal function). 

Theorem 5.2. To every constraint C n which is both a (primary) pre- and post-constraint at 
step n there is associated 

(i) a null vector of Hf- , 
(ii) a right null vector ofQ n , and 
(Hi) a left null vector ofQ n +i. 

Furthermore, C n generates a flow tangential to C n = C~ H which is a symmetry of the 
(effective) action. 

Proof. The equations (15. lip applied to the case ~C n = + C n give the associated right and left 
null vector to £l n and fi n +l respectively. These equations also show that the right and left null 
vector coincide with each other and is given by ^pr- Subtracting the right equation in (|5,9p 
from the right equation in (|5.10p . both applied to this case, yields 

f)(l n 

-Q-n^ = (5-13) 

which shows that this vector also defines a null vector of the Hessian (evaluated on the space of 
solutions) . 

Finally, {x l n ,C n } = ^tt, thus the flow is defined by a null vector of the Hessian and the 
Lagragian two-forms and hence generates transformations in flat directions at the extrema of 
the action. Given that C n is first class with all other constraints, it follows that the flow is 
tangential to the constraint surface C n = n C~ . □ 

Theorem 5.3. To every constraint C n := ~C n = + C n which is both a primary pre- and 
post- constraint at step n there is associated a gauge mode. 

Proof. From what has been said before, evolution from (n — 1) to n and from n to (n + 1) 
leads to an a priori and a posteriori free parameter A n = which is not fixed by any of the 
other pre- or post-constraints at n. We have to ensure, that this does not change if we evolve 
either further into the future or into the past, since pre- or post-constraints might propagate 
to n. However, all these evolutions could be described by effective actions, all of which, by the 
conjunction of theorems 13.11 and 14. 1 \ would lead to the same pre- and post-constraints (including 
the propagated ones). Thus, X n = fi n remains a free parameter throughout evolution because 
its corresponding constraint always remains first class. □ 
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In conclusion, although the constraints in the discrete do not generate the dynamics, cases 
(a)-(c) of section [5TT1 lead to some similarities (and dissimilarities) with the continuum situation: 

(a) Coinciding pre- and post-constraints are first class and generate symmetries of Hamilton's 
principal function. 

(b) Pre-constraints which are independent of the post-constraints and do not fix any A n (and 
vice versa with the fi n ) are first class, but do not generate any symmetries. 

(c) Pre-constraints which are independent of the post-constraints at n, yet which fix a priori 
free A n must be second class together with those post-constraints whose flows they fix. 

That is, a symmetry generating constraint is necessarily first class also in the discrete (case 
(a)). However, in contrast to the continuum, a first class constraint does not necessarily generate 
symmetries in the discrete (case (b)). We will better understand this case (b) in the context 
of propagating degrees of freedom in section 15.41 below. Finally, recall that in the continuum 
Lagrange multipliers can only get fixed if they are associated to second class constraints. Like- 
wise, a priori free parameters X n in the discrete can also only get fixed if they are associated to 
post-constraints which are rendered second class by pre-constraints (case (c)). 

5.4 Varying numbers of propagating degrees of freedom 

We shall now discuss role (iii) of the constraints, namely the classification of degrees of freedom 
into gauge modes and gauge invariant observables. 

In the continuum the number of physical degrees of freedom is tied to classifying constraints 
into first and second class. The (standard) discussion of the continuum corresponds to translation 
invariant systems, where, in particular, the number of physical observables is constant during 
time evolution. In contrast, in our general set-up, the constraint surface at a given time n 
depends on the initial and final time step, thus we also have to expect that the notion of the 
number of propagating degrees of freedom depends on this choice. 

5.4.1 Observables as propagating degrees of freedom 

The principal idea is to define observables by propagation of data. In the continuum (or trans- 
lation invariant systems) the classification of the constraints at a single instant of time into first 
and second class yields a complete characterization of propagating observable and gauge degrees 
of freedom [13]. On the other hand, in order to specify the meaning of propagating degrees of 
freedom in the discrete, we need two time steps; for the notion of propagation at the canonical 
level requires the global Hamiltonian time evolution map H n -\ '■ ^n-i ~~ ^ ■ 

As discussed in section 13.3.21 in the presence of constraints this map is a priori not well 
defined, rather, we need to specify a priori free parameters A n . Similarly, for the inverse evolution 

— > C~_ 1 a posteriori free parameters fi n -i must be specified. Varying these free parameters 
leads to orbits in the post- and pre- constraint hypersurfaces, respectively. The Hamiltonian 
time evolution map T-L n -i is well defined and invertible as a map from the space of orbits at 
(n — 1) to the space of orbits at n. In particular, these spaces are of equal dimension (see also 
theorem 13. ID . 

Consequently, these two spaces of orbits can be identified with each other via the Hamiltonian 
evolution map and specify the canonical data which can be either uniquely pre- or postdicted. 
These spaces, therefore, describe the propagating degrees of freedom of the evolution move 
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(re — 1) —7- n between the time steps (n — 1) and n. The number of propagating degrees of 
freedom will be defined as the dimensional of the space of orbits. 

Accordingly, we define pre-observables 0~_^ at step (n — 1) and post-observables 0£ at 
step n as functions which are well defined on the space of orbits, i.e. as fi n -i, A n -independent 
functions; via % n the pre-observables at (n — 1) uniquely predict the post-observables at n and, 
vice versa, the post-observables at n uniquely postdict the pre-observables at (n — 1). 

Theorem 5.4. The pre-observables 0~_ 1 Poisson commute (weakly) with the pre- constraints 
at (n — 1) and the post-observables 0+ Poisson commute (weakly) with the post-constraints at 
n, i.e. 

= 0, and {0,+ , + C r ra } = 0, VZ,r. 



Proof. The post-constraint + C™ generates a flow 

d+C? 



dp] 



d+C? d 2 S n d+C? d+p]d+C r r 



(n n + C n \ = — = J r (K -\A\ 

iP] ' r 1 94 c+ dxidxi dp? dxi d P f [ ■ > 

In (|5.1ip we showed that the gradient -gpf- ='■ (R n )i is a right null vector of the Lagrangian 

two-form. Such a null vector leads to non-uniqueness of x n in the image of the time evolution 
map, see section [3.3.21 This is infinitesimally described by x n + e\ r n (R n ) r where are the a 
priori free parameters. Thus, the first equation in (|5.14p shows that the constraint generates the 
corresponding flow for the configuration variables. The second equation shows that the post- 
constraint flow also reproduces the induced non-uniqueness for the post-momenta. Hence, the 
orbits discussed above are defined by the flow of the post-constraints. This flow is integrable, 
as the set of post-constraints is first class. 

A phase space function which commutes (weakly) with the post-constraints is constant along 
a given orbit and will not be affected by the non-uniqueness of the time evolution map. A similar 
argument holds for the pre-observables. □ 

The number of propagating degrees of freedom depends on initial and final step and the 
notion of observables is move dependent. This is illustrated by the following examples. 

Example 5.2. Let us reconsider examples 15.11 (A) - (C) of section [531 

(A) The evolution move — > 1, in which three field variables and three momenta evolve into 
the field variables and three momenta is regular. No constraints arise and the number of propa- 
gating degrees of freedom is six. The following move 1 — > 2 evolves three phase space pairs into 
two pairs. There is a pre-constraint at time 1 such that the number of propagating degrees of 
freedom for this move is four. For the effective evolution — > 2 a pre-constraint at arises and 
the number of propagating degrees of freedom coincides with the minimal number of the two 
moves, which is four. 

(B) In this case two phase space pairs at time evolve into three pairs at time 1 which subse- 
quently evolve back into two pairs at time 2. One finds a post- and pre-constraint at time 1 
and the number of propagating degrees of freedom for both moves is four. This also agrees with 



18 This is twice the number of propagating configuration degrees of freedom, which might be used as the definition 
of propagating degrees of freedom in the Lagrangian picture. 
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the number of propagating degrees of freedom for the effective evolution move — s> 2. 
(C) Three pairs at are evolving into two pairs at time 1 and then back into three pairs at 
time 2. We have a pre-constraint at and a post-constraint at 2 such that again the number 
of propagating degrees of freedom for both time steps is four which coincides with the number 
of propagating degrees of freedom in the effective evolution. 

The number N^f of independent propagating degrees of freedom between an initial step i 
and some final step /, is given by the dimension of the post-constraint hypersurface + Cf modulo 
the flow generated by the post-constraints. This agrees with the dimension of the pre-constraint 
hypersurface ~C{ modulo the flow generated by the pre-constraints. Hence, 

Ni-tf = 2Q — 2# (pre-constraints at i) = 2Q — 2# (post-constraints at /) (5.15) 

is the number of both the pre-observables at i and post-observables at /. This number coincides 
with the rank of the symplectic form restricted to the pre- and post-constraint surfaces discussed 
in theorem 13.11 However, is in general not the dimension of the reduced phase spaces at 

steps % or / because possible post-constraints at % and possible pre-constraints at / have thus 
far been ignored. 

5.4.2 The reduced phase space 

Let us define the reduced phase space at a given step n with i < n < f as the total constraint 
hypersurface C n := n C~ modulo any first class flows. As we shall see shortly, this reduced 
phase space will describe the degrees of freedom propagating from i through step n to /. 

Let us check how many of the Aj_> n propagating degrees of freedom of the move i — > n 
continue to propagate in the move n—¥f. We have to impose both post- and pre-constraints 
at n. There are three different cases that can arise, which we labelled by (a)-(c) in section [5TT1 

(a) Part of the (possibly re-organized) post-constraints coincide with part of the (possibly 
re-organized) pre-constraints. According to theorem 15. 2\ this part defines constraints 
that are first class with all other constraints at n and which, furthermore, generate gauge 
symmetries at n. Hence, for every such constraint the number of physical degrees of 
freedom decreases by two. However, this is automatically accounted for by the number 
Ni^ n of post-observables at n such that this case does not prevent any of the latter from 
continuing to propagate during n — > f. 

(b) We discussed in sections 15 . 1 1 and [5 . 31 also the possibility of a pre-constraint, that is indepen- 
dent of the post-constraints, yet is first class with all post-constraints (and automatically 
with all pre-constraints). Such a first class pre-constraint is related to an a posteriori 
free parameter jjL n on which the post-constraints do not impose any conditions. Hence, 
although this \x n is not a proper gauge mode and will have propagated during the move 
% — > n (as it is not a priori free), being a posteriori free, it will not continue to propa- 
gate during n —¥ f. That is, for any such first class pre-constraint, one canonical pair of 
the N^ n propagating modes of the move i —> n will cease to propagate for n — > /. (The 
analogous discussion holds for a post-constraint that is first class with all pre-constraints.) 

(c) Pre- and post constraints which are second class. For a second class post-constraint 
there is at least one pre-constraint with which this post-constraint has (on the constraint 
hypersurface) non-vanishing Poisson brackets. As the post-constraints generate the flow 
of the parameters A n and the pre-constraints generate the flow of the parameters fj, n , it 
follows that a corresponding pair of these variables (X n ,fJ-n) gets fixed by imposing the 
pre-constraint or the post-constraint, respectively. This can happen in two ways: 
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1 . The free parameters are determined by the right or left null vectors of the 

Lagrangian two-forms. Thus, despite having + C n ^ ~C n the null vectors may coincide 
and define the same free parameter X n = fi n . The latter is fixed to a certain value by the 
constraints, i.e. the a priori and the a posteriori orbits are intersected transversally by the 
constraint hypersurface. However, none of the Ni_> n post-observables at n is prevented 
from propagating further to /, because A n = [i n is both a priori and a posteriori free (and 
thus did not propagate in any direction in the first place). 

2. A n 7^ fj, n for the parameters associated to post- and pre-constraint, respectively. Again, 
both parameters are fixed by the constraints to certain values and do not decrease further 
the dimension of the reduced phase space. More precisely, the a posteriori free parameter 
/x n is only free for the move n — > f , however, is not a priori free for the move i — > n and 
thus propagates during the latter £J but does not continue to propagate during n —> f. The 
converse is true for the a priori free A n . Consequently, an observable pair, corresponding 
to fi n , propagates during i —> n and does not propagate further. However, during the 
fixing via the constraints, the data of this pair is transferred to another observable pair, 
corresponding to A n (which was not an observable pair for i — > n), that continues to 
propagate during n — >• f. Thus, none of the N^ n propagating data is 'lost' in this case. 

Let us choose a basis of constraints which separates the first and second class constraints. 
The respective numbers are defined independently of the choice of this basis (e.g., since the 
number of first class constraints agrees with the number of degenerate directions of the pull 
back of the symplectic form to the constraint hypersurface). 

Combining all of the above, the number of degrees of freedom that propagate from i via n 
to / is 

Ni^ n -).f = Ni^ n — 2 x ^t(pre-constraints of case (b) at n) 

= 2Q — 2 x #(lst class constraints at n) — ^(2nd class constraints at n) (5.16) 

which thus coincides with the dimension of the reduced phase space at n. Note that this 
dimension depends on the choice of initial and final time steps i and /. 

Theorem 15.41 has an immediate consequence for a special kind^l of observables: 

Corollary 5.2. Observables O n := 0+ = 0~ which are both pre-observables of the move 
n — > / and post-observables of the move i —> n at time step n, Poisson commute (weakly) with 
all constraints at time step n. Thus, the O n are elements of the space of functions which are 
well defined on the reduced phase space associated to time step n of the evolution i — > n — > f. 

Example 5.3. We consider the analogue of the 'no boundary' proposal [14j in Regge Calculus, 
as briefly discussed in [2] . Take an evolution move — >■ 1 from the empty triangulation at n = 
to the boundary of a single D-simplex, D > 3, at n = 1. The phase spaces To, V\ corresponding 
to a single simplex are totally constrained [2] and, hence, the corresponding reduced phase spaces 
are zero-dimensional. No information (i.e. 'lattice gravitons') propagates during —¥ 1 (a single 
simplex is flat and does not contain any 'gravitons') When performing further simplicial 
evolution moves [2] to a larger spherical hypersurface at step n (see figure [T4"|) and solving the 

19 That is to say, fj, n actually gets predicted by the initial data at step i which under Hi propagates to the 
post-constraint surface and thus is automatically compatible with the post-constraints at n. 

20 Notice that observables, in general, do not need to also Poisson commute with the second class constraints. 
This suggests to reconsider the interpretation of the 'graviton propagator' from spin foam models which has 
been derived from a single 4-simplex (e.g., see [37] )■ The semiclassical limit of spin foam models is expected to 
be dominated by Regge geometries. However, as argued, a classical Regge 4-simplex does not contain 'gravitons'. 
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intermediate equations of motion, one finds that V n is still totally constrained (this follows from 
theorems 13, II and 14. ip . Thus, the reduced phase space at n is also zero-dimensional. Indeed, no 
information can propagate from the empty triangulation at to any other spherical hypersurface 
ri, or, colloquially, "no information can propagate from 'nothing' to 'something'." 




Figure 14: Discrete analogue of the 'no boundary' proposal |14j . 



This does not imply that all spherical 4D Regge geometries are devoid of 'lattice gravitons' 
because the notions of observables and the reduced phase space are step dependent. Namely, 
consider, instead, a 4D evolution from a smaller (non-empty) intermediate spherical hypersur- 
face at some i > to a larger spherical hypersurface at some / <n (see figure [T5"]) . In general, 
the piece of 4D triangulation interpolating between these two hypersurfaces will contain curva- 
ture and will not lead to totally constrained phase spaces at i and /, such that propagation for 
i — > f is possible. 

In conclusion, the pre-constraints at i and the post-constraints at / of a given (effective) 
move i — > f determine what propagates from i to /, while the conjunction of pre- and post- 
constraints at the same step n determines what propagates through n. The classification of the 
constraints into first and second class only plays a role in determining the reduced phase space 
at a given step and the presence of gauge symmetries. 



5.4.3 Propagating degrees of freedom and local moves 

So far we have only considered observables as propagating degrees of freedom under global 
evolution moves. Let us now also consider the local evolution of a given time slice. This is 
described by momentum updating as given in section 14.21 

By locally evolving a given time slice forward and solving any equations of motion arising 
on the way, one always considers the propagation of information from some initial step ki onto 
the evolving slice (which, in general, depends on ki). That is, evolving a time slice forward is 
tantamount to studying the sequence of global moves ki — >■ k, ki — >• k + 1, ki —¥ k + 2, ... It is 
clear that the number of degrees of freedom that propagate from the fixed ki onto the evolving 
slice cannot increase, but at best remain constant and in general will decrease. This is the 
'physical' reason why the rank of the symplectic form restricted to the post-constraint surface 
of the evolving slice — which determines the number of propagating degrees of freedom — can at 
best remain constant, but will in general decrease (see theorem 14. ip . 



Figure 15: Schematic 
evolution from a smaller 
to a bigger spherical hy- 
persurface. 
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On the other hand, it is also possible to create new degrees of freedom at, say, k + 1 that only 
propagate from k + 1 onwards. These new variables can be viewed as new initial data which 
only become relevant at step k + 1. This highlights the fact that in general discrete systems 
with evolving phase spaces sufficient initial data for the entire evolution is only fully assigned 
in the course of evolution. 

We briefly discuss the specific role of the four types of local moves of section 14.21 

Type I & IV: These moves do not reduce the number of observables propagating from ki onto 
the evolving slice because they preserve the symplectic form (see theorem 14. ip . The type 
I move additionally introduces new a priori free variables at step k + 1. If these new data 
are not a posteriori free for some move k + 1 — >• k + x, they will propagate from k + 1 to 
k + x. 

Type II and III: These moves reduce the rank of the symplectic form (see theorem 14. ip . and 
thus the number of propagating observables on the evolving slice, by two for each of their 
pre-constraints that is independent of the post-constraints at k and is not rendered second 
class (case (b) of sections 15.11 15.31 and 15.4.2]) . 

Since some new 'initial data' may be a priori freely chosen at k + 1, non-uniqueness of 
solutions arises. The system is non-hyperbolic, because a fixed step k cannot predict/postdict 
the entire evolution; information propagating from it may eventually stop propagating without 
evolution in general breaking down. 

6 Conclusions 

The present work develops a general canonical formalism and constraint analysis for arbitrary 
variational discrete systems. In particular, it applies to discrete mechanical, lattice field theory 
and discrete gravity models with variational action principle. The formalism is equivalent to 
the covariant formulation, encompasses both global and local discrete time evolution moves and 
can handle both constant and evolving phase spaces. We have discussed the (non-) preservation 
of the constraints and symplectic structure in this framework and analysed the different roles 
constraints can assume in the discrete. In particular, constraints can 'propagate' on solutions 
such that the number of constraints at a fixed step n depends on the initial and final steps of 
evolution. We have shown that the sets of pre- and post-constraints at any step each form a 
first class Poisson-sub-algebra, and can only become second class when both sets are considered 
together. Gauge symmetry generating constraints are simultaneously pre- and post-constraints 
and thus necessarily first class. However, first class constraints can also arise which, in contrast 
to the continuum |12tll3|. do not generate symmetries but restrict the (physical) solution space. 
In analogy to the continuum, second class constraints lead to a fixing of free parameters. 

We have defined observables as propagating and gauge modes as non-propagating degrees 
of freedom. It is established that, in general variational discrete systems where the phase spaces 
vary from step to step, the meaning of a propagating degree of freedom depends crucially on 
the initial and final step between which one considers the time evolution. Nevertheless, the 
definition of pre- and post-observables is unambiguous and the general canonical formalism 
describing their propagation fully consistent. More precisely: 

• The notion of an observable at a single step n is in general not useful. 

• The number of propagating degrees of freedom of the evolution move i — > / depends 
on the initial and final time steps i and / and is determined by the pre-constraints at i and 
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the post-constraints at /: the N^f pre-observables at i Poisson commute (weakly) with 
the pre-constraints at i and the N^f post-observables at / Poisson commute (weakly) 
with the post-constraints at /. In general, for different i and / one has different numbers 
of propagating degrees of freedom. Thus, the number of observables generally varies. 

• The reduced phase space at a step n in an evolution i — > n — > f depends on i and / 
and coincides with the space of observables that propagate from i through n further to /. 
In particular, a totally constrained phase space at step n does not imply the absence of 
observables, but only the absence of observables propagating through n. 

• A meaningful counting of degrees of freedom is provided which, e.g., entails that the 
number of propagating degrees of freedom from a fixed initial step onto an evolving slice 
can only remain constant or decrease. 

For translation invariant systems, on the other hand, the determination of the space of observ- 
ables reduces to the well known procedure from the continuum. 

Systems with evolving phase spaces are generally non-hyperbolic, because a fixed step n 
cannot predict the entire future; information propagating from it (in both past and future) may 
eventually stop propagating without evolution in general breaking down. For instance, new 
'initial data' introduced during a local move k — > k + 1 can propagate from k + 1 onwards. Since 
some new 'initial data' may be a priori freely chosen at k + 1, non-uniqueness of solutions arises. 

In this article, we have so far worked with phase space extensions to suitably formulate 
systems with evolving phase spaces. It should be noted, however, that one may equally well 
define and work on a total phase space ~p tot : corresponding to all variables that ever become 
relevant in the entire discrete solution. For instance, in a Regge triangulation, V tot is the 
phase space associated to the entire space-time (and not hypersurface) triangulation of a given 
covariant solution. The dynamics on V tot proceeds then as follows: any variables that, at a 
given step re, are not relevant (either because they have not yet appeared in the evolution or 
have become 'bulk') are 'switched off' by the corresponding canonical constraints p n = and 
only the relevant variables at a given step are 'excited'. Their dynamics then simply proceeds 
by momentum updating. From this larger picture it is straightforward to return to the phase 
spaces V n that we used in the present article at step n: one only has to perform a phase space 
reduction on V tot using the first and second class constraints that appear at re. 

Aspects of the present general formalism allowed us in [2] to construct the canonical for- 
mulation of Regge Calculus, based on gluing single simplices at each time step onto the spatial 
hypersurfaces. This provides the attractive intuitive picture of the Regge dynamics as an evolu- 
tion of the triangulated hypersurfaces in a discrete 'multi-fingered' time through the full Regge 
solution, akin to the evolution of spatial hypersurfaces in canonical General Relativity. It also 
constitutes an algorithm to generate solutions, given a suitable set of initial data. 

As an open question for future research we would like to mention the issue of broken diffeo- 
morphism symmetry for discrete gravity [HI [18] . This implies that modes of the metric which 
are gauge in the continuum will actually be classified as propagating degrees of freedom in the 
discrete (if one considers a curved background solution). Thus, we cannot claim that all the 
propagating degrees of freedom (even on a regular lattice) necessarily map to the graviton (i.e. 
spin 2 ) modes of the continuum. The additional modes have been termed pseudo gauge modes 
[HUH], as one expects their dynamics to differ from the 'true graviton modes'. In particular, the 
correlations associated to these modes could be very short ranged compared to the correlations 
of the true graviton modes. The tools provided in this work allow to identify triangulated hy- 
persurfaces between which observables can propagate in principle. In the case of gravity these 
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observables will include both gravitons and pseudo gauge modes. As discussed, there are also 
hypersurfaces between which propagation cannot occur. It would be of great benefit to develop 
a systematic criterion to distinguish pseudo gauge modes from modes that are truly propagating 
also in the continuum limit. This certainly would improve the comparison of lattice results for 
the graviton propagator [37j to the continuum. Promisingly, the present formalism is amenable 
to coarse graining techniques under which pseudo gauge modes can be expected to average out. 

Finally, there are many interpretational questions regarding a possible quantization of this 
framework. Hilbert spaces associated to phase spaces of different dimension have to be mapped 
to each other under time evolution. In this context, the usual unitarity requirement for time 
evolution has to be replaced. In particular, the concept of dynamical cylindrical consistency 
|30| becomes essential, as it connects measures on Hilbert spaces based on coarser and finer 
discretizations. The 'general boundary formalism' [19|. I20j. on the other hand, offers useful tools 
for devising a meaningful probability interpretation for such evolving Hilbert spaces. Research 
in this direction is currently under way. 



A Relation between effective canonical and Lagrangian two- 
forms 



In this appendix we shall prove equation (|3.43j) of section 13.51 For simplicity and brevity, we 
shall assume that boundary data constraints are absent. We will only show this for Vt n under 
the post-Legendre transform as the proof for f2 n +i under the pre-Legendre transform proceeds 
identically. 

Prior to imposing the equations of motion the Lagrangian two-form f)3.5[) on Q n _i x Q n 
reads 



a, 



Differentiation of 



d 2 S„ 



jdx^_i A dx e n 



d 2 S n 



dx e n _ x dxl 

9S n (Xn_i, xi,x^ 



dx e n _ x dx\ 



~dx n 2 A dx n . 



(A.1) 



dxi 



yields 



d 2 s n d 2 s n dxi 



dx^dx 1 ^ 



dx^dx^ dx!^ 
which can be employed in 



0. 



d 2 S, n 



+ 



0. 



d 2 s n dxi 



dx t n dx^ l _ l dx^dx^ dx^_^ 



(A.2) 



dx!^ 



I dXyy "J" 



dx 



dx 



n-l- 



n-1 



By the conjunction of the above and assuming d t ^" t , is invertiblel 22 ! one obtains from (|A. 1[) 

ox n ox n 



the effective Lagrangian two-form on Q n ~i X Q n 

( d 2 S n d 2 S n ( d 2 S 



a, 



d 2 S n 



^2 

dxidx e n _ x dx^dx^ ydx^dx^ J dx\ t dx n _^ 

d 2 Sr, 



dx n ^ A dXy 



d x i-\d x n 



dx e n _ x A dx e n . 



(A.3) 



If this matrix is degenerate, one can factor out the degenerate directions and invert the resulting matrix, 
which, for simplicity, we shall not worry about here. 
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On the other hand, the symplectic form on T*Q n reads 

oj n = dx e n A dp n e + dx l n A dp™. (A.4) 

Using the embedding i n : (x%,p%) h> = Pli x n = XniPt = °) of T *Qn into T*Q n |l 

we can pull back the symplectic form (|A.4|) to obtain, 

& n := i* n U) n = dx e n A , (A.5) 

the symplectic form on the 'reduced phase space' V n := T*Q n . (Note that additional constraints 
C(Xn,p™,) = may occur on V n .) Direct computation shows that on-shell one, indeed, finds 

n n = (¥+S n )*Cb n . (A.6) 

Analogously, one shows that fJ n +i = (F~S n )*u! n . 



B Momentum updating as a canonical transformation 

In section HT2l we mentioned that it is also possible to formulate momentum updating alternatively 
as a canonical transformation. Let us briefly elaborate on this. Since momentum updating applies 
to all canonical pairs on the extended phase space equally, we can extend the corresponding time 
evolution map beyond the pre- or post-constraint surfaces defined by p k = or p k+l = to 
the full extended phase space without specifying the momenta at k beforehand. That is, the 
extended time evolution map H k defined on the full phase space at k reads for all variable pairs 
i = b, e, n, o 



n k+i _ k , 9S k+1 
x k+i — x k> Pi —Pi+jr^ ■ \ a - l > 

ax k+\ 



The gauge fixing conditions x£ = and x° k = x^ +1 (which were a free choice in the previous 
formulation of momentum updating) are thus built in from the outset. The image of (jB.ip is 
given by the full extended phase space at time k + 1. The extended momentum updating (jB.ip . 
which involves a generating function of the second kind 

G(x k ,p k+1 ) = Y / xl kPt 1 -Sk+i, 

i 

preserves the canonical two-form (of the extended phase space). Namely, by symmetry and 
antisymmetry and using Uk+i = dx\ +1 A dp^ +1 , one immediately verifies 

(H k )*LJ k+1 = dx\ A d (p k i + ^j 1 ) = dx\ A dp\ + d f k a + \, dxi A dxl = dx l k A dp\ = to k . 
\ ox k J ax k ax k 

Subsequently, one must impose the pre- and post-constraints p\ = and p k+1 = 0, respec- 
tively, by hand. Note that, as a consequence of x k = x k+1 , this leads to the post-constraints 

+ C k+1 = p k+1 — (which are the image of the pre-constraints p k = under (IB.ip ) and the 

pre-constraints ~C k = p k Q + 9S d k Ji 1 (which are the pre-image of p k+l = 0) on the extended phase 
space. However, the latter pre- and post-constraints need not in general be constraints on the 



23 Viewing the x^-\ in this case as 'external' parameters. 
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unextended phase spaces due to the general dependence of the action on both x k+1 anc ^ x °k- ^ n 
this case, the pre-constraints = and ~C k and, likewise, the post-constraints p k+1 = and 
+ C k+1 , in fact, are generally second class. Those ~C k and + C k+1 which are only constraints 
on the extended, but not on the unextended phase space (originating, after all, in the general 
condition x l k = x l k+1 ) can be viewed as gauge fixing conditions such that not only x k = x k+1 
and x° k = x k+1 , but also x k+1 and x° k are no longer free parameters. 

This implies preservation under time evolution of the induced symplectic forms on the com- 
bined gauge fixing and constraint hypersurfaces. On these gauge fixing and constraint surfaces, 
the extended momentum updating is equivalent to the previous version of momentum updating 
which was only formulated as a pre-symplectic transformation. 



C Proof of theorem 14.11 of section 14.21 

Proof. Consider type I. Take the symplectic forms at k and k + 1 on Vk^Vk+i 

u k = dx\ A dpi + dx\ A dp k e + dx k A dp k n , u k+1 = dx b k+l A dp k+1 + dx e k+1 A dp k e +l + dx k+1 A dp k+1 . 

Notice that p!^ = is both a pre- and post-constraint at k. Restrict the symplectic forms to 
the partial post-constraint surfaces defined only by the K p\ = at k and only the K post- 
constraints + C k+1 on the right in (|4.6p at k + 1 and denote them by K- k ,K k+1 . K, k+l is the 
image of K k under rjfc and vice versa. At time k employ the canonical embedding jk : JC k > Vk 

Jk ■ {x\,p k ,X%,p k 1 X k ) h-> (Xf.jpfr, X%,p k , XfrjP^ = 0). 

The pull-back of the symplectic form to K, k coincides with the symplectic form of the unextended 
phase space Vk-, 

(jk)*0Ok = dx k A dp\ + dx% A dp k . (C.l) 
For time step + consider the embedding j^+i : K k+1 <^4 Vk+i given by 

/ i) k-\-\ e fa I 1 n \ / b fa I 1 e fa I 1 n n \ 

Jk+l '■ (X k+1 ,p b ,X k+1 ,p e ,X k+1 ) (X k+1 ,p b ,X k+1 ,p e , X k+l ,p k+l = — — j. 

UX k+l 

The corresponding pull-back of the symplectic form to JC k +i rea( is 

(jk+iTuk+i = dx b k+1 A dp k+1 + dx e k+1 A dp k+l + dx n k+1 A f Sk+1 e dx% +l 

OX k+l OX k+l 

= dx\ +1 A dp k+1 + dx% +1 A d [p k + l - ||±i \ . (C.2) 

Using the coordinate form (j4.4H4.6p of t)k to pull back the pre-symplectic form ()C.2p at time 
k + 1 , we find 

t)Ujk+i)*Uk+i = dx b k Ad P k + dx e k Adp k = (j k )*u k . 
Any post-constraint that further reduces the rank of (|C.2|) must be of the formal 



+ r k+l I e b fc+1 dSk+l fc+1 

u I ^fc+D^fc+l'-Pe ~'^b 



0. 



24 Here and in the sequel an equation such as p k — + p k represents a shorthand notation for equation (|3.13p 
p k = + p k = + 9Sk ^ k ' X0 ^ for the post-momenta, where Sk is the (possibly effective) action for the (global) move 
— > k for some (here irrelevant) initial step 0. The analogous shorthand notation is used for pre-momenta. 



44 



By (|4.4H4.6|) . each of these pulls back one-to-one to a post-constraint at k, 



h * + ^k+l _ +r<k( e b k k\ 
■Jfc u — u {x k ,x k ,p e ,p b ) 

Hence, i) k preserves the total post-constraint surfaces, f)& : C k -+ C k+1 and (i k )*co k = hk( L k+i)*^k+l) 
where i k : C k <— >■ 7\ and t^+i : C^+i <— >■ 'Pfc+i are the embeddings of the post-constraint surfaces 
at and k + 1, respectively, into the extended phase space. 

Next, consider type II. The proof proceeds analogously by replacing the label n by o, but 
taking into account the presence of non-trivial pre-constraints ~C k . Denote by JC k the partial 
pre-constraint surface at step k defined only by the K pre-constraints ~C k in (|4.9|) . No other 
pre-constraints are created in the move. Denote by /CjJ", 1 the partial (post-)constraint surface 
at step k defined only by the K (post-)constraints p k+1 = which are the image of the ~C k 
under \) k . Let the corresponding embeddings be j k ■ JC k ^ Vk and j k+ i 

■ -+ Vk+i- In 

analogy to (fC31ICT2l) . one obtains 

(j k Tu k = dx\ A dp k b + dxt Ad(p k e + 9Sk+1 ^! (C3) 
(j k+1 )*u k+1 = dx b k+1 A dp k b +l + dx% +1 A dp k e +1 , (C.4) 

and, consequently, by @3H2J|, htOk+iY^k+i = {jkY^k- 

Any additional post-constraint at k + 1 that further reduces the rank of must 

be of the form + C k+1 (x b k+v x e k+1 ,p k+1 ,p k+1 ) ^ ^ ^ = 0. Using lj£7HP]> . it is pulled back 
one-to-one to a post-constraint at k 



r k ( T h <r e n k n k I dfffc+lQEfc i x %) | _ h * + r fc+lc 6 e fc+1 , 

u ^fc'^fc'Pb'Pe + j. , t — "Jfc ° l x fe+l; x k+vPb ' " 



p K =+p k 



0. 



which likewise reduces the rank of (jk)*^k- Hence, f)fc preserves all independent rank reducing 
constraints. Thus, we must have t) k : C k D /C^T — >■ CjJ, 1 and (i k )*ui k = M* k {^ k +i)*^k+i with the 
corresponding embeddings i k : C k n /C^ ^+ Pfc and ifc+i : C^" +1 ■— > "Pfc+i, where CjT'^fc+l are the 
total post-constraint surfaces at k and A; + 1, respectively. 

Consider type III. Again, p\ = constitute both pre- and post-constraints at k. Denote 
by Kt the partial post-constraint surface at k defined by the K p\ = only and by j k its 
embedding in P k . Similarly, denote by /CjJ" +1 the partial post-constraint surface at k + 1 defined 
by the K p k+1 = only and by j k+ \ its embedding in V k+ \. This yields 

(j k )*UJ k = dx\ A dp\ + dx e k A dp k + dx k A dp„ , 
(j fc+ i)*o; fc+ i = dz£ + i A + dx e k+1 A + dxl +l A , 

which coincides with the canonical forms of the unextended phase spaces V k ,V k +i. Furthermore, 
using (l4~T0ti4TT3i one finds 

J> a , , e A . / k , dSk+i \ , j__ n a j i d<Sfc+i 1 



h£( Jfe+1 )*^ +1 = dx b k A dpg + da:£ A d [p k e + + dx n k+1 A d , (C.5) 

= dx| A dp£ + dsg A dp k e +dx° k Ad = 0fc)*Wfc- 
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<9 2 S 

If the rank of ~ oa k t, 1 is K — k, there will arise k post-constraints 

OX k aX k+l 



C k+1 (x e k+1 ,xl +l ,p k+1 ) _+ „fc+i_££fc+i _ ( C - 6 ) 



from the second equation in (|4.12p which further reduce the rank of (jk+i)*0Jk+i by 2k. The 
pullback of ([01]) under d4TITOfll51) yields the identities t)* k + C k+l = + C k+l {x e k , a£ +1 , §§|±i) = 0, 

i.e. trivializes the constraints at fc and thus does not reduce the rank of (jk)*^k- However, in 
complete analogy the K pre-constraints 



C k (x k , x%,p k ) 



s k+1 = 



which arise from the second equation in (|4.13p and define /C fc reduce the rank of {jk)*^k by 2k 
and are promoted to trivial identities at k + 1. 



Any additional post-constraint + C k+1 (x b k+1 , x k+1 , p k+1 , p k+1 , p k+1 ) fc+ = at 

A; + 1 that further reduces the rank of (jk+i)* '^k+l is pulled back, employing (|4.10fJ3~T3j) . to a 
non-trivial post-constraint + C k (x b k , xt, x k+1 ,p k ,p k + d ^t 1 , gfn +1 ) = at fc which, by 

fc fe+l pk = -\-pk 

virtue of ()C.5p . likewise reduces the rank of (jk)*^k- in complete analogy, one verifies that 
also the reverse holds true such that fj^ defines a one-to-one mapping between the additional 
post-constraints at k and k + 1. The conjunction of the above implies \)* k {ik+i)*0Jk+i = (ikY^k, 
where tk '■ C k n JC k Pk and Lk+i '■ C k+l > Vk+i are the corresponding embeddings. 

Finally, the momentum updating map t) k °f type IV, given by ()4. 14t I4.15p . preserves the 
symplectic structure because no new pre- or post-constraints arise in the move and any already 
existing post-constraints at k are trivially mapped one-to-one to post-constraints at k + 1. □ 
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